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The  inviscid  and  viscous  stability  problems  are  addressed  for  a  boundary  layer  which  can 
support  both  Gortler  and  Crossflow  vortices.  The  change  in  structure  of  Gortler  vortices 
is  found  when  the  parameter  representing  the  degree  of  three-dimensionality  of  the  basic 
boundary  layer  flow  under  consideration  is  increased.  It  is  shown  that  Crossflow  vortices 
emerge  naturally  as  this  parameter  is  increased  and  ultimately  become  the  only  possible 
vortex  instability  of  the  flow.  It  is  shown  conclusively  that  at  suflSciently  large  values  of 
the  crossflow  there  are  no  unstable  Gortler  vortices  present  in  a  boundary  layer  which,  in 
the  zero  crossflow  case,  is  centrifugally  unstable.  The  results  suggest  that  in  many  practical 
applications  Gortler  vortices  cannot  be  a  cause  of  transition  because  they  are  destroyed  by 
the  3-D  nature  of  the  basic  state.  In  swept  wing  flows  the  Gortler  mechanism  is  probably  not 
present  for  typical  angles  of  sweep  of  about  20  degrees.  Some  discussion  of  the  receptivity 
problem  for  vortex  instabilities  in  weakly  3-D  boundary  layers  is  given;  it  is  shown  that 
inviscid  modes  have  a  coupling  coefficient  marginally  smaller  than  those  of  the  fastest  growing 
viscous  modes  discussed  recently  by  Denier,  Hall  and  Seddougui  (1990).  However  the  fact 
that  the  growth  rates  of  the  inviscid  modes  are  the  largest  in  most  situations  means  that 
they  are  probably  the  most  hkely  source  of  transition. 


‘Research  was  supported  by  the  National  Aeronautics  and  Space  Adirunistration  under  NASA  Contract 
No.  NASl-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 


§1.  Introduction. 

Recently  there  has  been  interest  in  the  effect  of  boundary  layer  growth  on  instability 
mechanisms  to  which  the  boundary  layers  are  susceptible.  Here,  we  shall  concentrate  on 
the  Gortler  vortex  mechanism  which  has  been  shown  to  occur  in  both  two-  and  three- 
dimensional  boundary  layer  flows  over  concave  walls.  Much  of  the  early  theoretical  work 
concerned  with  Gortler  vortices  addressed  the  problem  of  the  linear  stability  of  external 
two-dimensional  flows  over  such  concave  walls.  Early  contributions  were  made  by,  among 
others,  Gortler  (1940),  Smith  (1955)  and  Hammerlin  (1956).  Later  Hall  (1982a, 6)  argued 
that  much  of  this  early  work  was  fundamentally  flawed  for  all  the  analyses  invoked  the 
parallel  flow  approximation  (which  essentially  assumes  that  the  basic  flow  in  which  the  vor¬ 
tices  lie  is  independent  of  the  streamwise  co-ordinate  and  so  neglects  the  effect  of  boundary 
layer  growth).  This  approximation  enables  the  linear  stability  equations  to  be  expressed 
as  ordinary  differential  equations  but  Hall  illustrated  that  this  assumption  is  unjustifiable 
except  in  the  limit  of  small  vortex  wavelength,  and  indeed  this  is  the  explanation  for  the 
considerable  inconsistencies  in  the  results  of  the  previous  studies.  Additionally,  in  the 
case  of  small  vortex  wavelength  the  Gortler  instability  may  be  described  by  an  asymptotic 
structure  which  takes  account  of  boundary  layer'growth  in  a  rational  manner  and  thence 
the  parallel  flow  assumption  is  rendered  superfluous  in  the  only  situation  in  which  it  has 
any  relevance  whatsoever. 

Hall  (1982a)  examined  boundary  layer  flow  over  the  cylinder  y  =  0,  — oo  <  z  <  oo 
where  the  z-axis  is  a  generator  of  the  cylinder  and  y  measures  the  distance  normal  to  the 
surface.  The  x-  co-ordinate  measures  distance  along  the  curved  surface,  which  is  taken  to 
have  variable  curvature  (1/6)  k  (y)  where  6  and  I  are  lengthscales.  The  Reynolds  number 
Re,  the  curvature  parameter  6  and  the  Gortler  number  G  are  defined  by 

Uni  I 

Re=-}L,  6  =  -,  G  =  2Re^6,  (1.1) 

1/  b 

where  Uq  is  a  typical  flow  velocity  in  the  x  direction  and  u  the  kinematic  viscosity  of 
the  fluid.  Hall  (1982a)  investigated  the  flow  characteristics  when  Re  is  large  and  6  is 
small  such  that  in  the  limit  6  — >  0  G  is  held  fixed  at  an  0(1)  value.  By  scaling  the 
span  wise  co-ordinate  z  on  the  boundary  layer  thickness,  it  was  demonstrated  that  for 
a  non-dimensional  vortex  wavenumber  1  linearised  vortex  modes  within  a  two- 

dimensional  basic  flow  are  neutrally  stable  at  a  Gortler  number  G  =  go€~^  +  . . .,  where  go 
is  a  known  0(1)  constant  whose  precise  value  is  dependent  upon  the  properties  of  the  basic 
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flow  under  consideration.  For  larger  wavelengths  the  problem  is  fully  nonparallel  and  the 
linear  stability  equations,  which  now  take  the  form  of  a  set  of  partial  differential  equations, 
have  to  be  solved  numerically,  see  Hall  (1983).  This  paper  showed  two  significant  features 
of  this  nonparallel  flow  problem,  namely  that  the  ideas  of  a  unique  stability  curve  and 
of  unique  growth  rates  at  a  specified  downstream  location  are  inapplicable  to  the  Gortler 
problem  because  the  location  where  a  vortex  commences  to  grow  is  dependent  upon  the 
position  and  the  shape  of  the  the  imposed  disturbance. 

Questions  concerning  the  development  of  nonlinear  nonparallel  vortices  within  growing 
boundary  layers  were  addressed  by  Hall  (1988).  This  numerical  investigation  showed  that 
as  the  nonlinear  disturbance  evolves  the  perturbation  energy  becomes  concentrated  in  the 
fundamental  and  mean  flow  correction;  a  conclusion  consistent  with  the  weakly  nonlinear 
theory  of  Hall  (19826)  valid  for  small  wavelength  vortices.  It  is  well  known  that  Gortler 
vortices  set  up  in  an  experiment  conserve  their  wavelength  as  they  move  downstream. 
Since  the  boundary  layer  itself  thickens  it  follows  that  the  local  nondimensional  vortex 
wavenumber  becomes  large  as  the  vortex  develops.  Thus  the  small  wavelength  limit  in 
the  external  Gortler  problem  is  appropriate  to  the  ultimate  development  of  any  fixed 
wavelength  vortex  and  hence  sufficiently  far  downstream  in  many  flows  the  asymptotic 
work  of  Hall  (1982a,  6)  becomes  applicable. 

As  with  all  weakly  nonlinear  investigations,  the  results  of  Hall  (19826)  are  valid  only 
within  a  neighbourhood  of  the  point  where  the  imposed  perturbation  is  neutrally  stable. 
For  vortices  of  wavenumber  ^  3>  1,  their  development  downstream  of  the  point  of  neutral 
stability  is  governed  by  the  solution  of  a  pair  of  coupled  nonlinear  partial  differential  equa¬ 
tions  which  adopt  a  simple  asymptotic  structure  at  large  values  of  X,  where  eX  (X  =  0(1)) 
denotes  the  distance  of  the  vortex  downstream  of  the  neutral  point.  Formally,  for  large 
X,  Hall  Sc  Lakin  (1988)  showed  that  this  asymptotic  structure  could  be  used  to  deduce 
the  flow  configuration  for  fully  nonlinear  Gortler  vortices,  at  which  point  the  mean  flow 
correction  generated  by  the  presence  of  the  vortices  is  as  large  of  the  basic  (undisturbed) 
flow  itself.  The  vortex  structure  derived  by  Hall  &  Lakin  (1988)  essentially  consists  of  a 
core  region  in  which  the  vortex  is  concentrated  and  which  is  bounded  by  two  thin  layers 
in  which  the  vortex  activity  is  reduced  to  zero  exponentially.  Further  work  by  Hall  &: 
Seddougui  (1989)  has  shown  that  these  thin  layers  are  susceptible  to  secondary  instabili¬ 
ties  which  take  the  form  of  travelhng  waves  confined  within  these  layers.  This  theoretical 
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investigation,  together  with  the  weakly  nonlinear  account  of  this  form  of  secondary  insta¬ 
bility  described  by  Seddougui  &  Bassom  (1990),  provides  good  qualitive  agreement  with 
several  experimental  observations,  notably  those  of  Peerhossaini  L  Wesfreid  (1988a,  6). 

All  the  work  described  above  has  addressed  problems  which  arise  when  Gortler  vortices 
occur  in  two-dimensional  boundary  layers  but  in  many  practical  situations  in  which  Gortler 
vortices  are  known  to  arise  the  basic  boundary  layer  is  three-dimensional.  For  example, 
in  the  case  of  a  boundary  layer  flow  over  a  three-dimensional  obstacle  or  the  flow  over  a 
turbine  blade  the  three-dimensionality  of  the  basic  flow  is  potentially  crucial  and  should 
not  be  neglected.  Most  signiflcantly,  the  development  of  laminar  flow  airfoils  has  given 
rise  to  designs  which  have  two  areas  of  concave  curvature  on  the  lower  side  of  the  airfoil 
and  when  the  wing  is  swept  the  boundary  layer  flow  is  fully  three-dimensional  and  the 
previously  mentioned  analyses  are  largely  inapplicable. 

The  first  attempt  to  describe  this  three-dimensionality  effect  was  given  by  Hall  (1985) 
who  examined  the  Gortler  mechanism  in  flow  over  an  infinitely  long  swept  cylinder.  The 
results  obtained  were  quite  general  and  did  not  require  a  precise  description  of  the  particu¬ 
lar  boundary  layer  under  investigation.  Hall  (1985)  showed  that  it  is  the  relative  size  of  the 
crossflow  and  chordwise  flow  over  the  cylinder  wlxich  is  critical  in  determining  the  vortex 
structure.  He  demonstrated  that  as  this  ratio,  say  A,  was  increased  from  zero  the  first 
significant  change  in  the  vortex  structure  from  that  in  the  two-dimensional  case  occurs 
when  A  ~  iZe“  j,  where  Re  is  the  (large)  Reynolds  number  defined  in  (1.1).  Then  the  vor¬ 
tices  become  time  dependent  and  unlike  the  two-dimensional  case,  the  high  wavenumber 
modes  no  longer  have  vortex  boundaries  aligned  with  the  flow  direction.  Indeed,  as  the 
crossflow  increases  further,  the  neutral  vortices  have  axes  perpendicular  to  the  vortex  lines 
of  the  basic  flow.  The  neutral  Gortler  number  for  the  vortices  was  predicted  by  a  large 
wavenumber  asymptotic  analysis,  the  results  of  which  suggested  that  for  0(1)  values  of 
the  ratio  of  the  crossflow  and  chordwise  velocity  fields  the  Gortler  mechanism  is  ^Aobably 
unimportant  compared  with  Tollmien-Schlichting  and  crossflow  type  instabih.ies. 

An  investigation  into  the  effect  of  crossflow  on  weakly  nonlinear  Gortler  vortices  was 
made  by  Bassom  (1989).  Rather  than  concentrate  on  vortices  within  a  growing  bound¬ 
ary  layer  with  the  inherent  difficulties  of  providing  an  adequate  description  of  boundary 
layer  growth,  that  work  considered  vortices  within  a  curved  channel.  This  has  the  ad¬ 
vantage  that  the  problem  can  be  described  by  solution  of  o^'i'inary  differential  equations 
as  opposed  to  partial  differential  equations.  Bassom  (1989)  showed  that  for  this  problem 
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equilibrium,  streamwise  independent  weakly  nonlinear  vortices  could  only  persist  below  a 
certain  threshold  value  of  the  crossflow. 

One  aim  in  the  present  paper  is  to  extend  some  recent  work  by  Denier,  Hall  &  Sed- 
dougui  (1990)  (hereafter  referred  to  as  DHS)  into  a  three-dimensional  setting.  The  prin¬ 
cipal  aim  of  DHS  was  to  investigate  the  problem  of  providing  a  rational  theory  for  the 
receptivity  of  Gortln  vortices  and  specifically  how  the  vortices  may  be  triggered  by  wall 
roughness  elements.  In  the  course  of  this  work,  DHS  reconsidered  the  stability  of  a  vor¬ 
tex  in  a  high  Gortler  number  flow  by  implementing  a  linear,  spatial  stability  analysis.  In 
particular,  for  G  1,  they  considered  the  structure  of  vortices  in  two  different  wavenum¬ 
ber  regimes:  they  reworked  the  analysis  of  Hall  (1982a)  for  high  o(g*)  wavenumbers 
and  also  examined  0(1)  wavenumber  vortices.  This  latter  disturbance  mode  is  governed 
by  inviscid  equations.  By  considering  the  region  between  these  two  wavenumber  regimes, 
DHS  identified  a  new  structure  which  is  relevant  for  vortices  of  wavenumber  o(g»)  and 

which  has  the  property  that  the  vortices  are  trapped  in  a  thin  layer  of  thickness  o(g“  ») 
at  the  wall.  This  is  in  contrast  to  the  structure  identified  previously  by  Hall  (1982a)  who 
demonstrated  that  vortices  of  wavenumber  o(g*  )  are  confined  to  a  thin  layer  embedded 
within  the  boundary  layer  at  some  well-defined  position. 

Of  paramount  importance,  DHS  showed  that  within  the  o(g»)  wavenumber  regime 
there  exists  a  unique  most  unstable  Gortler  vortex  according  to  linear  stability  theory. 
Hence  the  most  dangerous  mode  within  a  two-dimensional  boundary  layer  has  been  found 
and  our  objective  now  is  to  extend  the  work  of  DHS  to  consider  the  effect  of  introducing 
crossflow  into  the  boundary  layer  flow. 

As  previously  mentioned,  DHS  conducted  a  purely  spatial  analysis  of  the  instability 
modes.  Here,  we  firstly  add  a  crossflow  into  the  situation  considered  by  DHS  and  assess  its 
effect;  we  continue  to  use  the  spatial  approach.  For  the  0(1)  wavenumber  modes  at  G  1 
we  show  that  crossflow  first  has  a  significant  effect  on  the  two-dimensional  results  once 
this  parameter  becomes  0(iZe~  jG^).  By  examining  the  properties  of  the  vortices  when 
subject  to  this  size  of  crossflow  we  can  identify  several  limits  of  note.  We  show  that  as  the 
crossflow  increases  the  vortex  structure  takes  on  an  identity  which  is  essentially  that  of  a 
crossflow  instability,  a  mechanism  first  investigated  by  Gregory,  Stuart  &  Walker  (1955). 
Additionally,  if  both  the  crossflow  and  the  wavenumber  of  this  inviscid  vortex  increase 
appropriately  we  can  obtain  a  mode  whose  structure  is  dominated  by  viscous  effects.  It  is 
found  that  when  this  occurs  the  ratio  of  crossflow  to  chordwise  flow  becomes  0{Re~  a  G» ) 
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and  the  (large)  wavenumber  of  the  vortices  is  that  considered  by  DHS  and  which  pertains 
to  the  most  unstable  Gortler  vortex  in  a  two-dimensional  flow.  Consequently,  we  have 
the  basis  for  explaining  how  the  stability  characteristics  of  the  modes  found  by  DHS  are 
changed  by  the  presence  of  crossflow.  We  demonstrate  that  the  introduction  of  crossflow 
into  the  problem  can  have  a  stabilising  effect,  at  least  according  to  our  linear  theory.  In 
particular,  whereas  in  DHS  it  was  shown  that  stationary  vortices  are  necessarily  unstable 
at  o(g*)  wavenumbers,  this  is  no  longer  true  once  crossflow  terms  are  introduced.  In 
addition,  for  certain  crossflow  values  there  exist  neutrally  stable  vortex  modes  whilst  at 
large  enough  values  of  the  crossflow  no  vortex  instabilities  induced  by  centrifugal  effects 
persist. 

We  then  develop  our  description  of  the  disturbance  motion  by  allowing  the  vortices 
to  be  time  dependent.  Unfortunately,  as  in  DHS,  for  all  the  problems  treated  here  the 
equations  which  determine  the  vortices  form  a  pair  of  coupled  ordinary  differential  equa¬ 
tions  with  associated  boundary  conditions  and  this  system  is  only  amenable  to  a  numerical 
solution  for  0(1)  parameter  values.  We  study  this  numerical  problem  for  a  number  of  par¬ 
ticular  parameter  choices  and  also  investigate  the  asymptotic  description  of  the  vortex 
structure  in  both  the  limits  of  large  and  of  small  vortex  wavelengths. 

The  procedure  for  the  remainder  of  the  paper  is  as  follows.  In  the  coming  section 
we  formulate  the  problem  at  hand  and  obtain  the  differential  equations  which  determine 
the  (essentially  inviscid)  structure  for  vortices  of  0(1)  wavenumber.  After  considering  the 
solutions  to  this  problem,  we  move  on  in  section  3  to  examine  the  viscous  structure  which 
describes  the  large  wavenumber  0(G» )  vortices  in  the  presence  of  crossflow.  At  this  stage 
we  include  both  the  crossflow  and  temporal  terms  described  in  the  preceding  paragraphs 
and  then  subsequently  we  study  particular  cases.  In  section  4  we  briefly  consider  the  zero 
frequency  crossflow  problem  and  in  the  ensuing  section  examine  the  effect  of  allowing  the 
disturbance  to  have  periodic  temporal  dependence.  We  then  proceed  to  investigate  various 
asymptotic  hmits  which  describe  the  behaviour  of  the  vortex  modes  in  both  low  and  high 
wavenumber  cases  and  conclude  with  some  discussion. 


§2.  Formulation. 

Our  aim  here  is  to  obtain  the  equations  w’hich  determine  the  structure  of  high  wavenum¬ 
ber  vortices  in  a  slightly  three-dimensional  boundary  layer.  As  in  Hall  (1985),  we  consider 


a  boundary  layer  flowing  over  the  cylinder  y  =  0,  — oo<2<ooso  that  the  2r-axis  is  a 
generator  of  the  cylinder  and  y  measures  the  distance  normal  to  the  surface.  The  x  co¬ 
ordinate  measures  distance  along  the  curved  surface,  which  is  supposed  to  have  variable 
curvature  (l/6)x  (f)  where  b  and  I  are  length  scales.  We  define  the  Reynolds  number  Re, 
the  Gortler  number  G  and  the  curvature  parameter  S  as  in  (1.1)  where  Z/q  is  a  typical  flow 
velocity  in  the  streamwise  direction.  The  Reynolds  number  is  assumed  to  be  large,  whilst 
6  is  sufficiently  small  so  that  as  S  — >  0  the  parameter  G  is  fixed  and  is  order  one.  Again 
following  Hall  (1985),  the  basic  three-dimensional  boundary  layer  is  taken  to  be  of  the 
form 

u=  Uo[u{X,Y),Re-^v{X,Y),Re'^X*w{X,Y)\  +  O  {re' ^))  ,  {2.1a) 


where 

X  =  j,  Y={^)Re-K  (2.16) 

and  A*  is  supposed  to  be  0(1).  To  compute  this  boundary  layer  profile  in  practice,  ii  and  d 
are  determined  by  numerically  integrating  the  two-dimensional  boundary  layer  equations 
and  the  crossflow  w  is  then  found  from  the  span  wise  momentum  equation. 

We  now  define  the  variable  Z  =  Re^z/l  and  perturb  the  basic  flow.  If  t  is  a  time 
variable  scaled  on  I/Uq  and  if  =  exp{iaZ)  then  we  disturb  (2.1a)  by  writing 


u  =  Uo[n  +  6U{t,X,Y)E,vRe-^  +  SRe' ^V{t,X,Y)E, 

X‘Re-^w  +  6Re-^W{t,X,Y)E]  (l  +  o(Re->))  , 


(2.2a) 


where  6  <C  1.  Similarly,  we  perturb  the  basic  pressure  p  by  writing 


p  =  p  +  Re  ^P{t,X,Y)E. 


(2.26) 


Substituting  the  expressions  (2.2)  into  the  governing  continuity  and  Navier-Stokes  equa¬ 


tions  and  linearising  yields  the  perturbation  equations 

Ux  YVy  iaW  =  0,  (2.3a) 

tZj  -|-  uUx  “I"  U Ux  "1“  vUy  “I"  V tiy  "I"  A  wiaU  =  Uyy  —  a  U,  (2.36) 

Vt  -t-  nVx  +  UVx  +  VVy  +  V Vy  +  X*  wiaV  -f-  Gx{X)uU  =  —Py  +  Vyy  —  a^V,  (2.3c) 

W,  4-  uWx  +  X'Uwx  +  vWy  +  X'Vwy  +  X*  iviaW  =  -iaP  -f  Wyy  -  a^W.  (2.3d) 
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The  scalings  (2.2)  are  those  appropriate  to  Gortler  vortices  (see  Hall  (1982a)  for 
father  discussion).  We  note  the  key  features  that  the  streamwise  velocity  disturbance  is 
asymptotically  larger  than  the  other  components  and  that  with  these  scalings  the  pressure 
term  drops  out  of  the  streamwise  momentum  equation  (2.36).  Further,  viscous  effects 
are  brought  into  the  perturbation  equations  (2.3)  by  the  choice  of  spanwise  and  normal 
scalings. 

At  this  stage  we  consider  the  spatial  inviscid  Gortler  modes  of  0(1)  wavenumber;  in 
the  context  of  the  viscous  equations  (2.3)  this  corresponds  to  taking  the  limit  G  3>  1.  DHS 
showed  that  for  A*  =  0,  a  =  0(1)  the  vortex  structure  is  found  by  considering  expansions 
of  the  forms 

U  =  {uq{x,Y)  +  G~  ^ ui{x,Y)  +  . .  exp  J  pdx'^  ,  (2.4) 

together  with  similar  expansions  for  G~^V,  G"»W  and  G~^P.  Substitution  of  (2.4)  in 
(2.3)  yields  the  problem  for  the  amplification  rate  (5 

VO  -  (2.5a) 

with  associated  boundary  conditions 

ao  =  0  at  y  =  0,  ?;o  — >  0  as  T  — >  oo.  (2.56) 

DHS  observed  that  the  most  unstable  spatial  mode  which  is  a  solution  of  (2.5)  may 
be  written  in  the  form 

VQ  =  uexp{-aY) ,  (2.5c) 

2 

which  is  valid  for  all  a.  Thus  we  have  the  somewhat  remarkable  result  that,  so  long  as  u 
vanishes  at  the  wall,  the  fastest  growing  inviscid  mode  has  growth  rate  independent  of  the 
basic  state. 

We  are  concerned  with  inclusion  of  sufficient  crossflow  in  the  problem  so  that  the 
governing  equations  (2.5)  are  altered  by  this  effect.  We  find  that  this  occurs  when  A*  = 
O(G^)  and  then  we  write 

{U,  V,  VF,  P)  =  {  Uo{Y)  +  ...,  GH-o(y )  +  ...,  Wo{Y)  +  ...,  GPo{Y)  +  . . .} 

X  exp  ^Gj  J  Pdx^  , 

(2.6a) 
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A*  =  Ig^. 


(2.66) 


At  this  stage  we  concentrate  on  a  purely  spatial  approach  so  that  on  substituting  (2.6) 
into  (2.3)  we  obtain,  in  turn, 


{pUo  +  iaWo)  +  ^  =  0, 

(/3u  +  iAatl))  Uq  +  Qy  Vq  =  0> 
{j3u  +  ilaw)  Vq  +  x^iUo  - 
(/3q  +  iXaw J  Wo  +  Xwy  Vq  =  —iaPo. 


(2.7a) 


(2.76) 


(2.7c) 

(2.7d) 


^From  (2.7)  we  can  easily  show  that 


(— Xaw'j  i-  Aau;)  +  \('U'yy^j  Vq  =  ii^XUUyVQ. 

(2.8a) 

The  above  equation  then  controls  the  inviscid  growth  of  vortices  in  a  weakly  three  - 
dimensional  boundary  layer;  it  is  of  course  the  appropriate  generalisation  due  to  wall 
curvature  of  the  well-known  Rayleigh  equation  for  unidirectional  flows. 

Equation  (2.8a)  needs  to  be  solved  subject  to  suitable  boundary  conditions.  Firstly, 
we  demand  the  typical  inviscid  condition  at  the  surface  of  the  cylinder  !/  =  0  so  that 


Vo  =  0 


at  y  =  0. 


(2.86) 


In  addition,  to  ensure  that  the  disturbance  is  confined  within  the  boundary  layer. 


Vo  — ^0 


as  Y  — *  oo. 


(2.8c) 


These  boundary  conditions,  when  added  to  (2.8a),  constitute  an  oigenproblem  for  the 
spatial  growth  rate  /?  in  terms  of  the  scaled  wavenumber  a  and  the  scaled  crossflow  %. 
We  notice  that  since  this  problem  is  a  localised  one  about  some  point  x  =  x,  we  may 
scale  the  curvature  parameter  x  out  of  the  probl<  m  by  appropriately  redefining  /?  and 
Consequently,  for  the  remainder  of  the  work  described  here  we  shall  take  x  =  1- 

To  fix  our  ideas  upon  some  definite  boundary  layer  flow  we  chose  to  consider  (2.8) 
when  u  is  given  by  the  Falkner-Skan  profile 

Q  =  /'(V):  /’"  +  //''  +  Ul-(/')^)  =0,  /(0)  =  /'(0)  =  0,  /'(oo):-l.  (2.9a) 
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This  flow  is  closely  related  to  the  steady  flow  in  the  boundary  layer  along  a  surface  of 
revolution  near  a  forward  stagnation  point,  see  Jones  &  Watson  (1963).  Cooke  (1950) 
showed  that  for  a  flow  over  an  infinite  yawed  wedge  at  zero  angle  of  attack  then  if  the 
streamwise  basic  flow  is  given  by  (2.9a)  then  the  spanwise  component  w  satisfies 

w  =  g:  g"  +  fg'  =  g(0)  =  0,  g{oo)  =  1.  (2.96) 

The  dependences  of  u  and  ®  on  Y  are  depicted  in  Figure  (2.1).  We  observe  that  the 
Blasius  profile  would  not  be  a  suitable  choice  for  u  and  w  for  then  we  retrieve  the  scenario 
described  in  Hall  (1985).  Hall  showed  that  in  flows  with  zero  pressure  gradient  and  with  Q 
and  w  therefore  linearly  related  the  whole  problem  of  Gortler  vortices  in  three-dimensional 
boundary  layers  becomes  degenerate  and  can  b^  reduced  to  a  two-dimensional  case.  Con¬ 
sequently,  we  needed  to  focus  our  attention  on  other  physically  realistic  three-dimensional 
flows  and  the  choice  (2.9)  was  motivated  by  these  requirements. 

Given  this  basic  flow  we  then  solved  (2.8)  for  a  variety  of  parameter  regimes.  We 
firstly  considered  the  effect  of  increasing  the  crossflow  A  on  the  vortex  growth  rate  Pr  {Pr 
denotes  the  real  part  of  the  complex  number  P).  .The  results  are  illustrated  in  Figure  (2.2) 
for  a  selection  of  A  values  and  we  can  make  some  immediate  deductions.  When  A  =  0  we 
of  course  return  to  the  case  of  DHS  and  P  is  given  by  (2.5c).  However,  as  A  is  increased 
then  for  most  prescribed  wavenumbers  a  the  corresponding  growth  rate  decreases  so  that 
the  effect  of  increasing  the  crossflow  is  to  stabilise  the  vortex  flow. 

For  small  crossflows  X  (less  than  some  critical  value,  say  A^  «  4.7)  we  observe  that 
although  the  vortex  is  more  stable  than  in  the  two-dimensional  base  flow  case,  we  still 
have  >  0  so  that  the  vortex  remains  unstable  for  all  a.  Once  X  reaches  the  value  X^ 
however.  Figure  (2.2)  indicates  that  there  exists  a  region  of  wavenumber  space  in  which 
vortices  can  no  longer  persist;  indeed,  for  X  >  X^  there  exist  cut-off  values  ai(A)  and  a2(A) 
such  that  the  vortices  only  exist  for  a  <  ai  or  a  >  02-  As  A  increases  yet  further,  we 
see  that  ai  decreases  and  02  increases  and  in  addition  the  growth  rates  pr  in  the  lower 
wavenumber  regime  grow.  In  Figure  (2.3)  we  have  illustrated  the  dependence  of  /?,  upon 
X  for  those  values  of  a  less  than  ai(X).  We  can  infer  from  this  figure  that  as  A  — »  00  the 
value  ai(A)  converges  to  some  non-zero  value,  say  aioo  ,  <md  this  behaviour  is  explained  in 
the  asymptotic  work  elucidated  shortly. 

Figure  (2.4)  shows  the  behaviour  of  the  imaginary  part  of /j,  /?j,  as  the  wavenumber  a 
and  the  crossflow  X  vary.  We  immediately  deduce  from  Figure  (2.4a)  that,  for  0  <  X  <  X^, 
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0i  is  negative  for  all  a  ar  ’  becomes  increasingly  more  negative  as  a  — >  oo.  In  addition,  we 
know  from  the  solution  {^.oc)  of  DHS  that  -=  0  for  all  a  when  we  consider  a  purely  two- 
dimensional  flow  (A  =  0).  For  values  of  A  >  A^,  see  Figure  (2.46),  we  find,  as  previously 
mentioned,  that  the  unstable  vortex  cannot  persist  in  some  wavenumber  range  which 
increases  in  extent  as  A  grows.  In  Figure  (2.46)  this  feature  is  indicated  by  the  gap  in 
the  possible  values  of  for  A  >  A^.  As  A  is  increased  yet  further,  see  Figure  (2.4c),  /3j 
becomes  progressively  more  negative  for  any  prescribed  wavenumber  a  and  here  we  have 
plotted  the  dependence  pi(n)  for  a  <  <ii(A). 

We  turn  now  to  consider  various  asymptotic  solutions  of  (2.8).  Firstly,  we  study  the 
case  A  ::§>  1,  a  <  ai;  i.e.  the  left-hand  solution  branches  in  Figures  (2.26),  (2.46)  and  those 
branches  sketched  in  Figures  (2.3),  {2Ac).  If  we  write 

l3  =  pl  +  o{l),  (2.10a) 


and  substitute  this  expansion  into  (2.8a)  we  obtain  at  leading  order 

-  2  V 

(-i/3a  +  aw)  (4^  -  -  {-iPu  -r  a®  )  Vq  =  0. 


(2.106) 


Since  curvature  efects  are  now  negligible  at  zeroth  order  this  is  a  standard  Rayleigh  equa¬ 
tion  and  we  now  have  a  connection  between  our  three-dimensional  vortex  flow  and  the 
work  of  Gregory,  Stuart  &:  Walker  (1955)  (hereafter  referred  to  as  GSW)  who  were  con¬ 
cerned  with  a  description  of  stationary  crossflow  vortices  formed  in  the  boundary  layer 
above  a  rotating  disc.  This  instability  is  due  to  the  inflectional  character  of  the  basic 
velocity  orofile  and  GSW  showed  with  china-clay  techniques  that  the  crossflow  instability 
took  the  fv>rm  of  a  regularly  spaced  pattern  of  spiral  vortices  which  was  stationary  relative 
to  the  disc.  Stuart  (in  GSW\  using  inviscid  theory,  suggested  that  the  instability  could  be 
associated  with  a  particular  inflectional  profile  in  which  the  inflection  point  coincided  with 
a  point  of  zero  velocity  somewhere  in  the  flow.  His  calculation  gave  the  predicted  number 
of  vortices  to  be  approximately  four  times  the  observed  \alue  of  about  30,  but  the  angle 
of  13°  between  the  axes  of  the  vortices  and  the  radius  vector  was  in  excellent  agreement 
with  their  experiments. 

We  can  appeal  to  Stuart’s  inviscid  analysis  and  note  that  there  is  a  solution  ot  (2.106) 
with  boundary  conditions  (2.86,  c)  for  which  a  is  real  and  P  purely  imaginary  if  there  exists 
a  point  y  ^  0  at  which 

(—zpu(Y)  +  aw(Y))  =  (  —  ipu  (Y)  +  aw  ^Y))  =  0.  (2.11) 
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Some  numerical  work  shows  that  for  the  basic  profiles  (2.9)  a  suitable  solution  of  (2.11) 
may  be  found  and  then  A/a  «  —0.85.  Consequently,  we  have,  recalling  (2.10),  neutrally 
stable  crossfiow  vortex  type  solutions  of  (2.8)  if  A  ^  1,  a  =  0(1)  and 


^  -0.85. 
aX 


(2.12) 


On  figure  (2.4c)  we  indicate,  for  each  crossflow  parameter  considered,  those  points  for  which 
(2.12)  is  satisfied.  We  can  conclude  that  the  asymptotic  result  (2.12)  is  very  accurate  for 
surprisingly  modest  values  of  A.  Once  conditions  (2.11)  &  (2.12)  are  satisified,  the  Rayleigh 
equation  (2.106)  has  regular  singular  points  at  Y  =  0  and  Y  =  Y  and  the  solution  of 
this  equation  with  boundary  conditions  (2.86)  forms  an  eigenvalue  problem  for  the  GSW 
wavenumber  a.  We  recall  that  earlier  we  observed  that  for  A  >  A^  the  vortex  mode  can  only 
persist  for  wavenumbers  a  less  than  ai(A)  and  greater  than  a2(A)  and  that  as  A  — >  oo,  ai 
tends  to  a  non-zero  constant  aioo .  It  is  now  easy  to  see  that  aioo  is  the  GSW  wavenumber 
a.  We  have  chosen  not  to  compute  the  precise  value  of  a  here  for  its  value  is  dependent  on 
the  choice  of  the  basic  boundary  layer  flow  and  our  present  purpose  is  to  concentrate  on  the 
underlying  description  of  vortices  in  three-dimensional  flows  as  opposed  to  concentrating 
on  profiles  (2.9)  in  particular.  However,  from  Figure  (2.3),  we  can  deduce  that  for  the 
basic  flow  (2.9),  a  «  1.3. 

Thus  within  the  context  of  in  viscid  stability  theory  we  now  see  the  relationship  between 
Gortler  and  crossflow  vortices.  For  a  two-dimensional  boundary  layer  which  is  centrifugally 
unstable  the  growth  rate  of  the  most  dangerous  mode  is  given  by  (2.5c),  the  exact  solution 
found  by  DHS.  When  the  crossflow  is  increased  from  zero  a  maximum  value  of  the  growth 
rate  develops  at  some  wavenumber  and  the  growth  rate  curve  to  the  right  of  this  point 
eventually  crosses  the  zero  growth  rate  axis.  Thus  at  a  finite  value  of  the  crossflow  a 
finite  band  of  unstable  wavenumbers  persists  to  the  right  of  a  =  0.  At  this  stage  we  can 
identify  this  unstal  s  band  of  modes  as  crossflow  vortices  since  at  even  larger  crossflows 
they  reduce  to  the  modes  of  GSW.  At  larger  values  of  the  wavenumber  an  unstable  band 
of  modes  persists  up  to  a  =  oo  for  all  values  of  the  crossflow;  it  is  perhaps  appropriate 
to  think  of  these  modes  as  the  remnants  of  the  Gortler  mechanism.  However  we  shall  see 
below  that  at  high  wavenumbers  viscous  effects  become  important  and  that  when  these 
effects  are  allowed  for  these  remnants  vanish  completely  at  large  enough  crossflows. 

Having  examined  the  asymptotic  description  of  modes  for  which  A>l,a  =  0(l)we 
next  consider  vortex  states  at  high  wavenumbers  a.  In  Figure  (2.5)  we  show  the  dependence 
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of  the  eigenvalue  0  of  (2.8a)  for  a  range  of  crossflows  in  the  cases  a  =  5  and  a  =  10.  We 
observe,  as  already  commented  upon,  that  increasing  A  tends  to  decrease  the  growth  rate 
Pr  and  we  expect  that  for  each  fixed  a  ^  1  there  exists  a  corresponding  crossflow  value 
at  which  the  vortex  is  completely  stabilised  by  the  crossflow.  This  deduction  is  consistent 
with  the  results  previously  illustrated  in  Figure  (2.2)  where  as  the  crossflow  increases  the 
wavenumber  range  over  which  vortices  are  unable  to  persist  steadily  widens.  We  can  also 
conclude  from  earlier  results  that  for  a  ^  1  and  A  not  large  the  vortex  flow  will  remain 
unstable  for  as  a  — >  cxd,  /?,  ~  as  when  A  =  0  and  it  is  not  expected  that  an  0(1)  increase 
in  the  scaled  crossflow  A  could  result  in  a  reduction  in  the  growth  rate  by  this  dramatic 
amount.  Consequently,  we  anticipate  that  near  neutral  vortices  will  be  described  for  a  ^  1, 

A  >  1. 

When  a  ^  1,  following  DHS  it  is  deduced  from  (2.8)  that  the  vortices  will  be  confined 
to  a  thin  region  of  thickness  0{a~  ^)  adjacent  to  the  cylinder  wall  y  =  0.  We  therefore  find 
it  convenient  to  write 

i  =  (lY,  (2.13a) 

and  suppose  that  around  y  =  0  the  basic  flow  velocities  have  regular  Taylor  expansions  of 
the  forms 

Q  = + ...,  w  =  y2ia~^(  +  +  . . . .  (2.136) 

When  A  is  gradually  increased  we  find  that  the  stability  characteristics  of  the  vortices 
are  first  altered  significantly  once  A  ~  0{a3).  If  we  write  A  =  a^A,  A  ==  0(1)  we  are  led 
to  the  expansion 

P  =  i  (a  =  /?o  +  +  . . .)  ,  (2.13c) 

with  Po  and  Pi  real.  If  in  turn  we  expand  Vq  —  Vqo  +  a'^Foi  +  •••  and  substitute  into 
(2.8a)  we  obtain  at  the  leading  two  orders  the  relationships 

Poyn  +  =  0,  (2.14a) 


and 


+  -  {ftoyi2  + 


)- 


+  -  {0om2  +  >^^22)  \pQy12  T  A/f22]  Voo  =  (mii)^^^')30, 


(2.146) 
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which  needs  to  be  solved  subject  to  Voo  =  0  at  ^  =  0  and  as  ^  — >  oo  (recall  that  we  scaled 
X  =  1).  This  problem  needs  a  fully  numerical  solution  but  here  we  merely  remark  that 
this  equation  has  a  critical  layer  structure  at  the  position 


2/3ipii 


1/U2  + 


(2.14c) 


and  that  this  structure  has  to  be  fully  analysed  in  order  that  the  numerical  scheme  be 
effected.  We  are  presently  considering  this  problem  in  further  work  and  hope  to  be  able 
to  report  on  this  in  the  near  future.  However  the  structure  outlined  above  is  certainly 
consistent  with  our  numerical  calculations  which  suggest  the  emergence  of  a  critical  layer 
structure  when  a  approaches  the  cut-off  wavenumber  ai. 

For  the  present  purposes,  we  can  infer  that  Gortler  vortex  states  require,  at  wavenum¬ 
bers  a  ^  1,  a  crossflow  of  size  0{a^)  in  order  to  have  a  significant  effect  on  their  stability. 
Further,  in  this  case  the  whole  vortex  structure  is  compressed  into  a  thin  0(a"^)  region 
close  to  the  cylinder  surface.  Of  course,  in  the  development  of  the  inviscid  equations  (2.7) 
which  describe  the  vortex  we  found  that  for  0(1)  wavenumber  vortices  in  a  crossflow  of 
size  0(0^)  (see  (2.66))  the  viscous  terms  are  negligible.  However,  as  a  — >  oo  the  extent 
of  the  structure  reduces  so  that  at  some  stage  we  expect  that  the  viscous  terms  become 
significant  and  no  longer  negligible.  Guided  by  the  asymptotic  results  (2.13)  &:  (2.14), 
study  of  the  fundamental  equations  (2.3)  suggest  that  viscous  terms  are  crucial  and  so  a 
new  regime  is  acheived  when 

a  =  0(G'«),  (2.15a) 

and  then 

A  =  o(g^)  and  p  =  O  {g^)  .  (2.156, c) 

Recalling  the  original  scalings  (2.4),  (2.6)  we  then  have  a  viscous  structure  in  operation  for 
0{G> )  wavenumber  vortices  in  crossflows  of  size  0(G» )  and  then  the  vortices  develop  on  a 
streamwise  lengthscale  of  0(G"  ^ ).  Interestingly,  DHS  identified  the  unique  most  unstable 
vortex  within  a  two-dimensional  base  flow  as  occuring  within  this  wavenumber  band  and 
developing  on  this  lengthscale. 

In  thi"  section  we  have  restricted  ourselves  to  calculations  pertinent  to  stationary 
vortices.  Much  of  the  above  analysis  can  be  modified  in  a  straightforward  manner  to  ac¬ 
count  for  time-dependence  of  the  vortices.  However,  our  objective  has  been  to  concentrate 
on  understanding  the  fundamental  connections  between  the  inviscid  vortex  modes  and  the 
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GSW  crossflow  instability  and  viscous,  high-wavenumber  vortex  structures.  There  are  nat¬ 
ural  non-stationary  counterparts  of  the  GSW  stationary  instability  (see  Bassom  &:  Gajjar 
(1988)  for  an  example)  and  these  will  match  with  time-dependent  inviscid  vortices  in  a 
manner  virtually  identical  to  that  described  here.  Further  remarks  relating  our  inviscid 
modes  to  the  crossflow  instability  mechanism  will  be  made  in  the  conclusions.  However,  we 
now  move  on  to  examine  the  salient  properties  of  the  high  wavenumber  viscous  structures 
suggested  by  scalings  (2.15). 
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§3.  The  viscous  modes. 


Following  the  discussion  of  the  solution  properties  of  the  principally  inviscid  modes 
considered  in  the  previous  section  we  now  examine  the  effect  of  crossflow  on  the  stability 
of  the  viscous  modes  whose  scalings  are  suggested  by  the  asymptotic  results  (2.15).  The 
wavenumber  of  these  vortices  is  predicted  to  be  O  (g»),  i.e.  is  that  found  by  DHS  for 
the  most  unstable  vortex  in  two-dimensional  flow.  For  no  crossflow,  DHS  showed  that 
the  vortex  is  then  confined  to  an  0(G“»)  thick  region  which  lies  immediately  adjacent 
to  the  cylinder  y  =  0  and  that  the  spatial  growth  rate  of  the  instability  is  0(G»)  in  the 
streamwise  direction.  We  have,  from  (2.156),  that  when  the  scaled  crossflow  A*  =  O  (g») 
(where  A*  is  defined  in  (2.2a))  it  is  to  be  expected  that  a  significant  deviation  in  the 
stability  characteristics  of  the  flow  from  those  which  occur  in  the  zero  crossflow  case. 
Then,  following  DHS,  we  write 

a  =  A:oG«  ,  (3.1a) 

where  A:o  is  0(1)  and  seek  disturbances  confined  to  the  layer  in  which  =  0(1)  where 

rp  =  IcqG’Y.  (3.16) 

Then  we  consider  modes  for  which  the  perturbation  velocities  (U,V,W)  and  pressure  P 
take  the  forms 

{U,V,W,P)  = 

(z7o  +  G'  « +  .  •  • ,  Vo  +  G' Vi  +  . . . ,  GhFo  +  G»  +  . . . ,  G»  Pq  +  G>Pi  +  ...) 

X  exp  j^G*  J  {poiX)  +  G-  hiiX)  +  ...)  dX  -  iG^  J  (Do(0  +  G"  TDi(t)  +  . . .)  dt 

(3.2a) 

and  we  fix  the  crossflow  parameter  at  0(G»)  so  that 

A*  =  AG^  (3.26) 

Here  the  unknowns  Uq,  Ui,  Vq,  . . . ,  Pq,  Pi,  . . .  are  functions  of  0  alone. 

We  substitute  (3.1)  and  (3.2)  into  (2.3)  and  compare  like  powers  of  the  large  parameter 
G.  In  this  process,  we  need  to  expand  the  basic  flow  about  T  =  0  and  so,  as  in  (2.13),  we 
assume  that 

n  =  ^lniX)Y  +  i/ii2(X)y2  +  lfii3iX)Y^  +...,  (3.3a) 

Z  D 
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W  =  ^2l{X)Y  +  ^H22{X)Y^  +  l^i23{X)Y^  +  . 

Z  D 

V  =  o(y2),  for  r  <  1. 

i,From  leading  orders  in  (2.3a)  we  have 

PoUq  +  ikoWo  =  0, 

PiUq  +  PiUo  +  ^‘0^^  +  iArolFl  =  0, 
dxp 

and  similar  manipulation  for  (2.36)  yields 

PqIJ'II 


L  kQ  fcr^ 


ko 

iX 


‘-O 


2A:qPii 


+  *^21  =  0, 


Uo  =  f^Vo, 

Kq 


(3.36) 

(3.3c) 

(3.4a) 

(3.46) 

(3.5a) 

(3.56) 


and 


Idrp'^  '  ^0  ^0  2fcQPii 


iX 

X 


ey 

(M11A‘22  -  M2IMI2) 
iX  3  Pifiup^  ,  P2^^n^  I 


rr  T/ 

^1  -  1X^1 

Kq 


TT  .  ^12^1/ 

^0  +  1:3^  ^0- 

Kq 


(3.5c) 


On  considering  the  spanwise  momentum  equation  {2.3d),  at  leading  orders  we  just 
retrieve  a  trivial  combination  of  (3.4a),  (3.5a)  and  (3.56).  The  third  order  terms  imply 
that 

d2 

^  k  f 


L  dip^ 


iUo 

-  1  +  -p, - — 


iX 


^0 


^0 


2/cq/xii 


{fJ'llfJ’22  -  /^21M12) 


-  -7:2- 


lA  ^  ,3  ,  ,  /?2MllV' 

^(/illM23-M21M13)V'  +  — TTi —  +  - - 

Z/Cn  Kn 


Lepiifcg 


k^i 


^0 


VFo  +  — pr-^0  + 

^0  ^0 


(3.6) 

Eliminating  17i  and  Wi  between  (3.46),  (3.5c)  and  (3.6)  yields  an  equation  for  the  leading 
pressure  term  Pq  of  the  form 

dVo  .  dVo  (  iX^  ,  .  ,  Pwu\ 

dip 


Po  =  ^o-rry  -  ko-rr  +  I  772 —  {MnM22  -  M21M12)  +  .2  j  ^0 

Vfcg/rii  fcg  / 


dip' 


illo  dVo  PlfJ'llip  dVo 
+  - - 771 


iX 


ko  dip 


fcg  dip  k^mi 


-  TT—  (/^ll^22  -  /^21M12) 


dVo 

dip 

(3.7a) 
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Finally,  leading  order  terms  in  the  remaining  momentum  equation  (2.3c)  give  that 


ifio  Pmiii’ 

^0  ^0  2A;q//ii 


(3.76) 


so  that  on  eliminating  Pq  between  (3.7a,  6)  and  simplifying  we  obtain 


iQo  fiifiiip 


^0  2A;Q^ill 


(/ni/^22  -  Ai21/il2)t^' 


+  73 - (/^11^22  -  P21/^12)  ^'^0  =  -  Uq. 

kof^n 

It  is  of  interest  here  to  note  that  if  we  had  chosen  the  crossflow  to  have  been  of  the 
more  obviously  significant  size  0(G» )  then  the  quadratic  terms  in  in  the  above  equations 
would  be  replaced  by  line'^r  ones.  Then  j3i  can  be  redefined  to  remove  these  terms  so  that 
such  a  weak  crossflow  does  not  have  a  significant  efect  on  the  viscous  modes.  To  summarise 
thus  far,  the  leading  order  vortex  quantities  Uq,  Vq  satisfy  the  sixth  order  system  (3.56), 
(3.8).  To  simplify  this  pair  we  introduce  the  parameters  k,  i3,  U  and  V  defined  by 

^0  =  (x^ll)  ’  =  Uq  =  huU  and  Vq  =  x‘ ■  (3.9a  -  d) 

Further,  if  we  write 


IIq  =  and  A 


A  inilfi22  -  ^J■12^J■2l) 

2^Sx'» 


(3.9e,/) 


and  substitute  (3.9)  in  (3.56)  and  (3.8),  we  obtain  our  final  equations  to  determine  U  and 
V,  i.e. 


iCl  iXip^ 

^3  ^3 


iX-ip^X  f  cP  A  22 A  -  xJ^O 

k3  ~  k3  J  \  dih^  ^  +  13  r-3  ’ 


(3.10a) 


(— - 


^  _  iCl  pip  iXip^ 


^  -  V 


(3.106) 


We  need  to  impose  suitable  boundary  conditions  on  the  solutions  of  (3.10).  First,  in 
order  to  ensure  that  the  disturbance  is  confined  to  the  thin  0(G~  •  )  zone  close  to  the  wall 
we  demand  that  the  perturbation  quantities  decay  exponentially  as  ip  — >  oo.  Additionally, 
to  satisfy  the  requirement  of  zero  disturbance  velocities  on  the  wall  xp  =  0,  U  =  F  =  0 
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on  ^  =  0  and  the  continuity  equation  (3.4b)  further  imphes  that  dV/drp  =  0  on  =  0- 
Hence,  in  all,  our  six  boundary  conditions  are 

=  r  =  —  =  U,  at  V’  =  0.  (3.11a) 

dij' 

dV 

U  =  V  = - y  0  as  0  - >  oo.  (3.116) 

dtl) 

Now  (3.10)  and  (3.11)  constitute  an  eigenproblem  to  be  solved  for  the  scaled  vortex  fre¬ 
quency  fi,  scaled  spatial  growth  rate  p  and  scaled  crossflow  A.  We  note  that  once  the 
eigenfunctions  U,  V  are  determined  we  can  retrieve  the  span  wise  velocity  component  of 
the  vortex  by  applying  (3.4a)  and  (3.9). 

The  equations  solved  by  DHS  are  precisely  (3.10)  with  fl  =  A  =  0.  They  demonstrated 
that  for  all  k  we  then  have  /3  >  0  so  that  the  corresponding  vortex  flow  is  unstable.  Further, 
as  k  — >  0,  p  ~  k^  and  as  k  — >  oo,  p  ~  The  behaviour  as  k  — >  0  matches  with 

(2.5c)  which  determines  the  growth  rate  of  0(1)  wavenumber  vortices  in  two-dimensional 
flow.  Since  p  — >  0  as  both  k  — >  0  and  as  k  — — »  oo  there  is  a  most  unstable  vortex  at 
some  intermediate  value  of  k  which  DHS  calculated  to  be  =  0.476  with  corresponding 
P  =  0.312.  The  dependence  of  P  upon  k  as  found  by  DHS  is  illustrated  in  Figure  (3.1). 
Our  interest  is  with  non-zero  D  and  A,  but  before  we  consider  any  specific  calculations 
we  make  one  observation  concerning  (3.10).  We  see  that  if  (fi,/3.  A)  is  an  eigenvalue  set  of 

(3.10)  (3.11)  then  so  is  (  — D*,/5’,— A)  for  real  crossflows  A  and  where  an  asterisk  denotes 
complex  conjugate.  Evidently  we  can  therefore  restrict  our  attention  to  postive  crossflow 
parameters  and  do  so  for  the  remainder  of  the  paper.  As  a  first  detailed  examination  of 

(3.10) ,  (3.11)  we  will  study  the  effect  of  introducing  the  crossflow  A  whilst  maintaining  a 
zero  frequency  D —  this  case  is  treated  in  the  forthcoming  section. 


§4.  Solution  of  (3.10)  for  the  zero  frequency  case. 

We  first  consider  the  solution  of  governing  equations  (3.10)  with  associated  boundary 
conditions  (3.11)  for  the  steady  vortex  problem  (D  =  0).  Naturally,  this  problem  had  to 
be  tackled  numerically  and  the  solution  technique  used  in  practice  was  based  on  a  method 


18 


suggested  by  Malik,  Chuang  ^  Hussaini  (1982).  In  the  zero  frequency  work  described  in 
this  section  we  solved 


( 


^  ^  f  _  \ 

( _  1  _ 

^  d^’2  k^  k^  ^ 


2iX  -  it-.'U 

(4.1a) 

!■'  + - V  =  -- — 

(4.16) 

by  prescribing  values  for  the  real  scaled  wavenumber  it  and  the  crossflow  A.  Then  (4.1) 
with  its  associated  boundary  conditions  (3.11)  w'ere  solved  for  the  complex-valued  spatial 
growth  rate  /3.  Recalling  the  definitions  (3.2a),  (3.9b)  we  can  see  that  when  Re(p)  >  0  the 
vortex  is  unstable  and,  conversely,  for  Re(^)  <  0  the  disturbance  is  stable.  As  remarked 
already  we  can,  without  loss  of  generality,  restrict  ourselves  to  considering  A  >  0. 

Mahk  et.  al.  (1982)  described  a  boundary  value  method  for  solving  eigenvalue  prob¬ 
lems.  In  these  methods  the  differential  equations  to  be  solved  are  reduced  to  a  set  of  linear 
algebraic  equations  using  either  a  finite  difference  discretisation  or  a  spectral  representation 
and  the  eigenvalues  are  found  by  solving  the  characteristic  determinant  of  a  generalised 
eigenvalue  problem.  The  particular  method  devised  by  Malik  et.  al.  uses  a  fourth  order 
accurate  (Euler-Maclaurin)  finite  difference  scheme  with  nodal  points  distributed  so  as 
to  resolve  any  singular  layers.  These  authors  used  their  scheme  to  examine  the  temporal 
and  spatial  stability  of  a  three-dimensional  compressible  boundary  layer  flow  over  a  swept 
wing.  For  further  detials  of  this  numerical  method  the  reader  is  referred  to  Malil  et.  al. 
(1982). 

To  implement  the  solution  procedure  for  (4.1)  we  decided  to  perform  a  sequence  of 
runs  with  a  gradually  increasing  crossflow  A  and  a  fixed  wavenumber  k  =  ^.  This  choice 
of  k  is  not  significant,  and  we  duly  performed  further  tests  at  other  values  of  k  in  order 
to  ensure  that  the  results  were  qualitively  similar  to  those  we  describe  here.  For  zero 
crossflow  we  were,  of  course,  just  repeating  the  calculation  of  DHS  presented  in  Figure 
(3.1)  for  which  0  is  real  and  the  vortex  is  unstable  for  all  k  >  0.  For  non-zero  crossflows 
the  spatial  dependence  parameter  /?  becomes  truly  complex  and  some  solutions  of  (4.1)  are 
shown  in  Figure  (4.1).  Here  we  have  fixed  k  =  ^  and  w’e  observe  that  as  A  increases  Re(0) 
decreases  so  that  the  vortex  becomes  progressively  stabilised.  When  A  =  0,  ^  =  0.313 
(from  DHS)  and  Re(0)  vanishes  when  A  ~  0.406,  at  which  point  lrn(0)  aj  —0.91  (see 
Figure  (4.16)).  As  the  crossflow  increases  yet  further  i?e(/5)  becomes  negative  and  hence 
we  have  the  situation  of  a  stable  vortex  flow.  The  form  of  behaviour  of  p  when  f  ^  is 
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typical  of  that  for  other  wavenumbers  and  so  we  may  conclude,  at  least  for  these  stationary 
vortices,  that  the  effect  of  increasing  the  crossflow  is  to  stabilise  the  flow. 

VVe  then  extended  our  computations  by  fixing  the  crossflow  parameter  A  and  solving 
(4.1)  for  a  selection  of  wavenumbers  k.  The  results  obtained  are  illustrated  in  representative 
Figures  (4.2a,  b).  Here  we  plot  the  dependence  of  Re{p)  upon  k.  In  Figure  (4.2a)  we  have 
fixed  A  =  0.414,  i.e.  a  httle  larger  than  the  value  at  which  the  vortex  of  wavenumber 
k  =  0.5  becomes  neutrally  stable.  It  is  clearly  seen  that  although  Re{P)  <  0  for  =  0.5, 
there  is  a  range  of  wavenumbers,  approximately  0.37  <  k  <  0.44,  over  which  Re{^)  is  still 
positive.  Consequently,  at  this  crossflow  value  there  remains  an  interval  of  wavelengths  for 
which  the  stationary  vortex  is  unstable.  However,  from  Figure  (4.25)  it  is  plain  that  when 
A  is  increased  further  (here  A  =  0.434)  we  have  Re{P)  <  0  for  all  wavenumbers  and  hence 
when  the  crossflow  becomes  this  large  the  flow  is  completely  stabilised.  More  detailed 
calculations  showed  that  this  occurs  whenever  the  crossflow  is  greater  than  approximately 
0.416. 

There  are  two  natural  asymptotic  regimes  to  investigate  within  the  context  of  the 
present  problem.  The  first  deals  with  the  large  wavenumber  limit  k  - — ►  oo.  We  do  not 
address  this  limit  here  because  it  also  arises  when  making  large  wavenumber  studies  of  the 
temporal  problem  dealt  with  in  section  5.  Consequently,  we  shall  focus  our  attention  on 
k  — >  oo  there.  Second  is  the  case  of  small  scaled  wavenumber  k  — >  0.  We  can  develop 
a  formal  asymptotic  solution  of  (4.1),  (3.11)  in  this  limit  but  rather  than  concentrate  on 
small  wavenumbers  for  stationary  vortices  alone  w^e  shall  consider  a  slightly  wider  class 
of  flows  and  shall  investigate  small  wavelength  vortices  with  small  frequencies.  Of  course 
we  can  recover  the  asymptotic  solution  for  stationary  vortices  without  formal  difficulty 
but  also,  much  more  importantly,  considering  small  wavenumber  vortices  of  the  type  just 
described  will  allow  us  to  tie  this  analysis  in  with  results  for  the  temporal  case  in  section 
5. 


§.^.f  Solutions  of  (3.10)  for  A:  <C  1. 

Guided  by  the  eigpnsolutions  of  the  problem  (4.1)  for  small  k  it  is  found  that  for 
fc  “C  1,  fi  =  0{k'^)  the  solution  vortex  of  (3.10)  is  neutrally  stable  for  a  scaled  crossflow  of 
size  0(k^).  To  demonstrate  this  we  suppose  that 

n  =  nk^  +  ...,  \  =  Ck^+...,  (5  =  iDk^  +  ....  (4.2) 
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Here  we  have  assumed  that  ^  is  purely  imaginary  so  that  D  is  real  and  with  the  scalings  as 
in  (4.2)  the  solution  structure  divides  into  a  number  of  distinct  regions.  For  the  moment 
we  take  fl  >  0  as  we  will  subsequently  find  that  a  modified  asymptotic  structure  is  required 
for  the  case  <  0. 

Close  to  0  =  0  we  write 

7lj  =  0e,  V  =  Vq  +  ...,  U  =^k-^Uo+  ....  (4.3) 


Substitution  of  (4.2)  &:  (4.3)  into  governing  equations  (3.10)  shows  that 


( 


-  in  -  iD6  -  ice 


0 


d^Vo 


+  2iCVQ  =  -eUa, 


/  de2 

-  ih  -  iDe  -  Uo  =  Vo- 


(4.4a) 


(4.4b) 


We  expect  the  disturbance  to  be  largely  confined  within  this  O(fc’)  layer  and  so 
we  demand  that  as  B  ■ — >  oo  the  solution  decays.  Inspection  of  (4.4)  then  imposes  the 
conditions 

1  1 

as  9  — >  CO.  (4.5a) 


V^cc-  +  ..., 


Uo<x-^  +  ... 


Additionally,  boundary  conditions  (3.11a)  demand  that 

dVo 


Uo  =  Vo  = 


dd 


0 


at  0  =  0. 


(4.5b) 


Now  v4.4),  (4.5)  together  form  an  eigenvalue  problem  for  the  0(1)  real  constants  C  and  D. 
Before  considering  this  further,  we  need  to  check  that  the  disturbance  decays  correctly  so 
as  to  ensure  that  the  boundary  conditions  (3.11b)  applicable  as  0  - — >  co  are  indeed  met. 

Inspection  of  (3.10)  subject  to  the  scahngs  (4.2)  suggests  that  when  0  =  0(1)  the 
velocity  components  of  the  vortex  behave  according  to 


U  =  k-U  +  ...,  V  =  k^v  + 

using  (4.3)  and  (4.5a).  Substitution  of  (4.6)  into  (3.10a)  yields  the  equations 


(4.6) 


and 


-iCrp^U  =  t>. 


(4.7a) 


(4.7b) 
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To  match  with  {4.5a)  requires  V  oc  ^  as  ^  — >  0  and  the  solution  of  (4.7a)  with  this 
property  is 

V  ^  (4.8) 

3 

where  E  is  some  constant  and  K^,{xp)  is  the  modified  Bessel  function  of  order  i/.  Solution 
(4.8)  has  the  required  exponential  decay  as  0  — >  oo  and  so,  therefore,  U  also  decays 
exponentially  (see  (4.76)).  Hence  we  have  obtained  the  complete  asymptotic  description 
of  the  solutions  of  (3.10)  (3.11)  for  neutral  modes  when  “C  1,  fi  =  0{k'^ ),  fl  >  0. 

Technically,  to  complete  our  description  of  this  limit  the  solution  of  (4.4),  (4.5)  is 
desired  for  various  values  of  fl.  We  have  not  done  this  here  for  this  aspect  is  not  within 
our  principal  objectives.  The  asymptotic  work  above  demonstrates  that  for  fc  <C  1  only 
a  small  scaled  crossflow,  Oik^)  is  required  to  completely  stabilise  a  vortex  of  positive 
frequency  0(^7),  This  is  an  exceptionally  tiny  frequency  and  so  in  practice  the  scaled 
frequency  parameter  (l  defined  in  (4.2)  is  likely  to  be  large.  It  can  be  shown  (and  details 
are  available  from  the  authors)  that  for  >•  1  we  have  C  —  0{Cl* )  and  D  =  0(fl‘)  and 
this  information  matches  with  further  asymptotic  limits  to  be  discussed  in  the  following 
section. 

The  above  describes  the  asymptotic  solution  of  (3.10)  (3.11)  for  ^  <C  1  and  small 

but  positive.  We  can  show  however  that  for  Q  negative  a  completely  different  structure 
must  come  into  play,  and  in  order  to  illustrate  this  revised  case  we  will  work  with  stationary 
vortices  (f2=0)  although  our  analysis  may  be  easily  modified  to  account  for  |fi|  •C  1,  <  0. 

In  the  stationary  case  we  need  to  solve  (4.1)  for  k  <C  1.  We  write 


P  —  ik^P  +  ...,  A  =  fc’A  +  ..., 

V  =  Vb-t-...,  U  =  k^Uo  +  ..., 

in  the  region  where  0  =  0(1)-  Substituting  these  expansions  in  (4.1)  yields 

-i  (/30  +  A0^)  -  l)  ^0  +  2iXVQ  -  -0f7o, 

(/f0  +  A02)  Uq  =  Vo, 

so  that  2  - 

(/30  +  A02)  ^  ^  -  Vo)  -  2A  (^0  +  A02)  Vo  -  0  Vo, 


(4.9) 


(4.10a) 

(4.106) 

(4.10c) 
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which  is  an  eigenproblem  for  the  quantities  3  and  A  when  we  impose  the  inviscid  condition 

Vq  - >  0  as  ip  - >  0  and  demand  that  Iq - >  0  as  p  - >  oo.  For  neutral  modes  we 

need  0  real  and  from  the  numerical  work  illustrated  in  Figures  (4.1)  k  (4.2)  we  saw  that 
for  positive  crossflows  A  in  (4.1)  the  neutrally  stable  voitices  satisfying  (4.1)  typically  had 
lm{0)  <  0.  Consequently,  in  terms  of  (4.9)  we  anticipate  that  the  eigenvalues  0,  A  of 
(4.10c)  are  such  that  0  <  0,  X  >  0  and  so  (4.10c)  has  a  critical  la3'er  structure  surrounding 
ip  =  —/3/A.  Of  course,  the  critical  layer  problem  (4.10c)  is  closely  related  to  (2.146)  (in  fact 
it  is  just  a  scaled  version  of  (2.146))  which  governs  the  large  wavenumber  inviscid  modes 
studied  in  section  2.  As  discussed  there,  in  order  that  a  numerical  solution  on  (4.10c)  be 
possible  a  detailed  analysis  of  the  critical  laj'cr  structure  is  required  and  which  we  will 
consider  in  a  future  article.  However,  the  pertinent  result  from  this  critical  layer  problem 
is  that  it  represents  the  description  of  how  the  neutrally  stable  vortex  modes  are  affected 
as  we  pass  from  the  parameter  regime  of  section  2  to  those  considered  in  this  section. 

Although  the  above  description  is  specific  to  stationary  vortices  it  is  easily  generalised 
to  account  for  disturbances  of  frequencies  fl  =  Oik^)  and  negative.  The  effect  of  the  non¬ 
zero  frequency  is  to  introduce  two  extra  constants  into  the  critical  laj'er  equation  (4.10c) 
but  the  overall  properties  of  the  problem  are  virtually  unaffected.  Again  inclusion  of  these 
non-zero  frequencies  provides  scope  for  further  discussion  of  the  critical  layer  problem. 

To  summarise  thus  far,  we  have  shown  that  for  small  frequencies  the  effect  of  increasing 
crossflow  in  the  fundamental  equations  (3.10)  with  boundary  conditions  (3.11)  is  to  tend 
to  stabilise  the  vortex  flow.  Interestingly,  we  have  also  demonstrated  that  for  small  non- 
dimensional  wavenumbers  the  structure  of  the  neutrally  stable  modes  is  dependent  largely 
on  the  sign  of  the  frequency  of  the  modes.  For  frequencies  fi  >  0  the  vortex  structure  is 
dictated  by  the  solution  of  eigenproblem  (4.4)  whereas  for  fl  <  0  we  obtain  a  critical  layer 
structure  which  matches  with  the  large  wavenumber  inviscid  modes  described  by  (2.14). 
Having  studied  the  low  frequency  problem  we  now  move  on  to  consider  larger  frequencies, 
indeed  we  examine  the  case  for  which  the  scaled  frequency  il  =  0(1). 
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§5.  The  vortex  problem  (3.10)  for  0(1)  frequencies. 

Here  we  examine  the  properties  of  Gortler  vortices  of  wavenumber  0(G»)  and  fre¬ 
quency  0(G»)  in  the  presence  of  crossflow  of  size  0(G»).  In  this  case,  of  course,  the 
scaled  parameters  Cl,  P  and  A  in  (3.10)  are  all  order  one  quantities.  For  the  purpose  of  this 
section  we  will  concentrate  almost  exclusively  on  neutrally  stable  states  for  the  results  of 
the  previous  sections  have  suggested  how  the  stability  of  the  flow  will  be  altered  by  chang¬ 
ing  particular  parameter  values  which  give  rise  to  neutral  nrodes.  To  solve  for  neutral 
modes  we  again  used  the  method  described  in  section  4  subject  to  a  few  slight  changes. 
We  then  specified  the  0(1)  frequency  Cl  and  the  wavenumber  k  in  (3.10)  and  treated  these 
coupled  ordinary  differential  equations  as  an  eigenproblem  for  the  crossflow  parameter  A 
and  the  spatial  development  quantity  /5  which  was  supposed  to  be  purely  imaginary.  In 
this  way  we  ensured  that  the  resulting  eigenvalues  corresponded  to  neutrally  stable  modes. 

Some  sample  results  are  illustrated  in  Figures  (5.1)  —  (5.3)  where  we  show  the  depen¬ 
dence  of  the  crossflow  A  and  the  parameter  on  the  scaled  wavenumber  k  for  the  three 
frequencies  fl  =  1,  0  and  —1.  Briefly,  we  can  see  that  neutral  modes  appear  to  be  possible 
over  the  complete  range  of  wavenumbers  k  and  that  for  large  k  the  crossflow  needed  to 
produce  neutrally  stable  modes  is  quite  small.  However,  a  striking  differerrce  between  the 
neutral  modes  is  observed  depending  on  whether  fl  >  0  or  D  <  0.  In  the  former  case  (see 
Figure  (5.1a)),  as  k  — >  0  A  tend  to  an  0(1)  value  whereas  in  the  latter  case  A  — >  0  as 
k  — >  0  (Figrrres  (5.2a,  5.3a)).  A  discussion  of  the  asymptotic  properties  of  these  two  cases 
is  given  later.  We  note  here  that  on  recalling  the  results  of  section  4  together  with  the 
zero  crossflow  calculation  of  DHS  we  can  infer  that  for  crossflows  greater  than  the  neutral 
mode  values  orr  Figures  (5.1a  —  5.3a)  the  corresponding  vortex  motions  are  stable  whilst 
the  opposite  is  true  for  lesser  crossflows. 

We  can  be  confident  that  the  results  shown  in  Figures  (5. 1-5.3)  are  sufficient  to  demon¬ 
strate  that  increasing  crossflow  again  stabilises  these  non-stationary  vortex  modes  in  a 
similar  way  as  was  described  for  the  stationary  case  in  section  4.  As  before,  there  exist 
several  regimes  in  which  asymptotic  description  of  (3.10)  k.  (3.11)  is  possible.  We  first 
concentrate  on  properties  of  this  eigenproblem  in  the  high  wavenumber  limit,  k  ^  1. 
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§5.i  The  high-wavenumber  (A:  ^  1)  limit 


In  this  limit  we  are  concerned  with  obtaining  an  asymptotic  solution  of  the  system 

(3.10) ,  (3.11)  for  k  1,  =  0(1).  We  can  initially  consider  the  problem  in  the  context 

of  neutrally  stable  modes  in  the  spirit  of  the  calculations  just  described.  It  has  been  seen 
that  for  low  wavenumber  modes  the  vortex  structure  appears  to  be  critically  dependent 
on  the  sign  of  the  frequency  f)  but  inspection  of  Figures  (5. 1-5.3)  suggest  that  this  is 
unlikely  to  be  the  case  for  high  wavenumber  modes  and  consequently  there  is  no  need  to 
specify  U  further  than  taking  it  to  be  of  order  unity.  An  idea  of  the  expected  structure 
can  be  gleaned  by  studying  Figure  (5.4).  Here  we  present  the  eigenfunctions  of  (3.10)  &: 

(3.11)  when  this  pair  of  equations  is  solved  for  neutral  modes  with  11  =  0  and  k  ^  1.35. 
Unfortunately,  for  values  of  k  much  larger  than  this  value  we  found  that  it  was  difficult 
to  obtain  satisfactory  nunrerical  conve..gence  but  the  evidence  of  Figure  (5.4)  and  close 
inspection  of  the  governing  equations  are  sufficient  in  order  to  strongly  suggest  that  the 
disturbance  is  confined  to  an  0(1)  thick  region  which  moves  away  from  ^  =  0  as  — >  oo. 
If  we  write 

li)  =  ck^  +‘y;,  (5.1a) 

p  =  tiPok^  +  ^hk-^+ ...)  ,  (5.16) 

and 

A  =  AqA'"  ^  +  All'"  ^  4- . . . ,  (5.1c) 

together  with  the  eigenfunction  expansions 


u  =  Uo  +  k~^Ui  +  V  =  kh'o  +  k-^Vi  + 


(5.1d) 


and  substitute  (5.1)  into  (3.10)  we  find  the  foUowing.  At  leading  orders  we  obtain 


00  =  — Aqc, 


(5.2) 


and,  at  next  order, 


-  1  +  ipQip  -  i  {Pic  +  Aic^) 


*■  dip'‘ 


-  1 


1  +  iporp  -  i  {Pic  +  Aic^) 


Uo  =  -cUo, 
Uo  =  Uo. 


(5.3) 
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By  redefining  ip  —  tp  —  {{Pic  +  Xic^)/Po)  we  simplify  these  equations  to 


The  precise  solutions  of  (5.4)  (5.5)  are  relatively  unimportant-  the  crucial  conclusion 

we  can  draw  from  this  asymptotic  work  is  that  for  an  0(1)  frequency  we  only  need  a  very 
small  0{k~  )  crossflow  to  stabilise  modes  with  wavenumber  k  ^  1.  In  passing,  we  notice 
that  the  eigenproblem  (5.4)  may  be  used  to  deduce  other  results.  In  particular,  we  have 
the  result  that  the  description  of  the  asymptotic  problem  k  - — y  oo  is  independent  of  the 
order  one  value  of  fl;-  a  result  which  is  in  keeping  with  the  numerical  evidence  presented 
earlier.  Indeed,  the  above  is  valid  for  aU  frequencies  H  =  ci(fc^)  so  that  in  particular  it  is 
valid  for  the  high  wavenumber  limit  of  the  zero  frequency  work  examined  in  section  4. 

We  now  make  some  comments  about  the  low  wavenumber  limit,  k  — ►  0. 


§5.2  The  low  wavenumber  (k  1^  limit  for  >  0 


We  have  already  observed  that  different  eigenproperties  of  (3.10),  (3.11)  are  expected 
for  fc  «C  1  and  that  the  sign  of  Cl  is  crucial  in  determining  these  properties.  Consequently, 
here  we  take  fl  >  0  and  examine  the  case  fl  <  0  later  (the  zero  frequency  case  has  already 
been  dealt  with).  When  fl  >  0,  fc  «C  1  the  eigenfunctions  of  (3.10)  are  concentrated  in  a 
thin  region  near  ip  —  0,  see  Figure  (5.5).  In  this  figure  we  illustrate  the  eigenfunctions  f7, 
V  for  the  case  fl  =  1,  =  1.3  x  10~^.  This  siiggests  seeking  a  multi-zoned  asymptotic 

solution  structure  of  (3.10)  for  k  <€.  1,  0.  =  0(1),  and  we  now  show  how  this  may  be 
achieved.  Initially,  we  suppose  that  the  crossflow  A  and  spatial  dependence  p  assume  the 
forms 

A  =  Aq  +  Aifc  *  +  X2k  ’  +  .  . . , 

P  =  i  [Pok^  +  Pik*  +  P2k  +  .  . .]  , 


2G 


(5.6a) 

(5.65) 


and  we  recall  that  we  are  seeking  neutral  disturbances  so  that  the  0(1)  constants  Aq,  Ai,  . . . , 
are  real.  After  experimentation,  we  can  infer  that  the  disturbance  is  largely 
trapped  in  an  0(A’*)  thick  region  which  is  located  at  a  distance  0{k^)  from  =  0- 
Guided  by  this,  when 

xj)  —  (5.7a) 

say,  with  x/xq  =  0(1),  we  need  to  expand 

U  =  Uok-  =  +  Uik-  ^  +  U2+  V  ^  Vo  +  Vik^  +  V2k^  +  ....  (5.7b) 


Substituting  (5.6),  (5.7)  in  (3.10)  yields  that  in  this  region  (refer  to  it  as  zone  I), 

^  +  /?oV'0  +  Aoti’o  — 


Po  +  2V’oAo  =  PiVo  +  =  0. 

If  we  define  Cl  —  /320O  +  A20o  ~  (Pil^^o)  write 


(5.8a) 
(5.8b,  c) 


we  deduce  that  within  zone  I  the  leading  order  disturbance  quantities  Uq,  K)  satisfy 


,  '.2] 

— T  —  xQ  —  i\otl>  I 


+  2iAoVb  =  —i’oUo, 


(5.9a) 


d2  -  ^2j 

— ^  -iQ  ~  iXo4’  ,  Uo  -  Vq. 


(5.96) 


For  the  disturbance  to  be  effectively  confined  to  zone  I  we  require  the  solutions  to 
decay  as  ^  — *  ±00.  We  can  show  that  the  potentially  most  dangerous  mode  is  even  in  xp, 
and  in  addition  it  is  easy  to  prove  that 


Vo  cx  i,  Uqcx  — - 


as  |  ^;|  - >  cx). 


(5.10) 


In  the  course  of  his  investigation  into  Gortler  vortices  at  0(G*)  wavenumbers  in 
three-dimensional  boundary  layers.  Hall  (1985)  considered  the  eigenproblem  (5.9),  (5.10). 
He  found  that  the  relevant  eigenvalues  are 


00 

—  =  4.71  and  -^  =  -3.59.  (5.11a,  6) 


27 


Even  though  we  now  have  enough  information  to  determine  Aq  and  (3q  via  {5.8a,  b)  and 
(5.11a),  we  must  check  that  this  zone  I  solution  obtained  above  can  be  satisfactorily 
matched  to  the  requisite  boundary  conditions  (3.11).  Since  zone  I  is  fixed  at  a  distance 
0{k^)  from  the  wall  ^  =  0  it  is  natural  to  define  a  second  region,  zone  II,  by 

xp  =  A->6>.  (5.12) 

By  (5.76)  and  (5.10),  within  II  the  eigenfunctions  must  expand  according  to 

V  =  VQk*  +  . . . ,  [7  =  +  •  •  •  ■  (5.13) 

Substitution  of  (5.6),  (5.12)  &  (5.13)  into  the  governing  equations  (3.10)  yields,  on  appli¬ 
cation  of  the  results  (5.8),  that 

(0  -  to?  ^  -  2vo  =  0,  (5.14a) 

and 

jAo(^  -  t6o)^wo  =  -VO-  (5.146) 

Strictly,  of  course,  we  need  to  consider  the  regions  6  <  ipo  a^nd  0  >  to  separately,  and  we 
shall  concentrate  on  the  former  first. 


The  region  d  <  to 


We  first  recall  that  the  boundary  conditions  (3.11a)  need  to  be  satisfied  at  ^  =  0 
and,  guided  by  result  (5.10),  w'e  solve  (5.14a)  subject  to  the  conditions  vq  — >  j  -f  . . .  as 
9  — >  to,  VQ  =  0  at  ^  =  0.  This  yields 


ti  -  {to  -  of 

to {to  -  0) 


(5.15a) 


and,  by  (5.146), 

«  [to  -  {to  -  6')^] 
uo  =  - r - -3- 

Ao0o  (''^0  -  d) 


(5.156) 


We  see  that  (5.15)  satisfies  only  two  of  ?io  =  vq  =  ^  =  0  at  0  =  0  and  to  force  the  third 
condition  to  be  met  we  now  add  a  viscous  layer,  zone  III,  at  ^  =  0.  If  in  this  layer  we 
define  the  0(1)  co-ordinate 


t  = 


(5.16) 
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and  the  expansions 


U  =  UQk*  +  . . . , 


V  =  VQk*  +  . . . , 


then  to  satisfy  (3.11a)  and  to  match  with  (5.15)  we  require 


(5.17) 


UQ  =  VO  =  — —  =  0  at  ^  =  0 

d(. 


(5.18a) 


3i  ;  3^ 

+  •  •  •  ,  VQ - y  +  ....  as  ^  - >  oo. 


n 

Substituting  (5.16),  (5.17)  into  (3.10)  yields 


^  =  1:572  - 'fij  “0  =  ”0. 


(5.185) 


(5.19) 


The  solutions  of  (5.19)  satisfying  (5.18)  are 


3(l-i)  f  ,  3i(l-0  3i^ 

L  2m0Q'\/2?T  -I  ^  2  /  ip^Clv^ 

where  =  zTi  and  H  is  chosen  so  that  ao(0)  =  0-  In  this  way  we  have  shown  that 
the  solution  proposed  so  far  satisfies  the  required  conditions  at  0  =  0.  We  now  need  to 
reconsider  the  0(^3 )  zone,  zone  II,  where  xj)  —  k^O  and  9  >  V'o-  Here  expansions  (5.13) 
and  equations  (5.14)  again  apply  but  this  time  we  need  to  solve  (5.14)  so  that  uq  — >  }  +  •  •  ■ 
as  9  — >  and  uo  — >  0  as  0  — >  oo.  The  required  solution  is  simply 


and  then 


{9  -  xpo)  ’ 


Ao  {9  -  i/'O)" 


(5.20a) 


To  complete  the  structure  we  are  led  to  consider  a  fourth  zone,  I\',  where 
and  from  (5.20)  this  means  that  the  eigenfunctions  have  behaviours  of  tin'  forms 


V  —  vQk  U  =  iiok  *  4-  . . . . 


(5.21) 
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►  0.  Substitution 

(5.22a) 

(5.226) 
0  together  with 

Do  oc  V’^A'i(V’),  (5.23) 

where  K^i'ip)  is  the  modified  Bessel  function  of  order  u.  Clearly,  iiQ  also  has  exponential 
decay  as  ip  — »  oo  and  so  we  have  a  complete  asymptotic  solution  of  (3.10),  (3.11)  for 
a  =0(1),  k  <C  1  and  p  imaginary.  The  full  four  zoned  structure  of  this  solution  is 
summarised  in  Figure  (5.6).  Returning  to  (5.8a, 6)  and  (5.11)  yields  that 

Ao^«0.46n«,  V'O  ^  l-47n»  and  /?o  ^  •  (3.24) 

Hence  we  conclude  that  for  small  wavenumber  vortices,  a  crossflow  of  approximately 
0.46Q*  +  ...  is  required  to  stabilise  the  flow.  We  note  also  that  as  — >  0  in  (5.24) 
our  predictions  match  with  the  solution  of  the  <C  1,  =  0{k^)  problem  described  in 

section  4.1.  Having  now  obtained  the  solution  structure  for  k  1,  U  >  0  we  turn  now  to 
consider  the  second  case,  fi  <  0. 

§5.5  The  low  wavenumber  (k  1)  limit  for  (7  <  0 

Having  deduced  the  difference  in  behaviours  of  solution  properties  of  (3.10),  (3.11) 
when  fc  <C  1  depending  on  the  sign  of  17  we  now  study  the  case  17  <  0.  From  Figure 
(5.2)  we  conclude  that  the  most  fundamental  difference  between  the  two  cases  is  that  for 
17  <  0  both  p  and  A  tend  to  zero  as  k  — >  0  iu  order  that  neutrally  stable  modes  may  be 
obtained.  In  order  to  study  the  structure  here  we  show  the  eigenfunctions  of  the  problem 
(3.10)  when  12  =  —1  and  k  =  8.7  x  10“^,  .see  Figure  (5.7).  Clearly  the  eigenfunctions 
are  confined  to  a  thin  region  centred  away  from  the  wall  0  =  0  and  suggests  that  some 
critical  layer  type  structure  is  operational.  This  i.;  in  agreement  with  the  low  frequency 


To  match  with  (5.20)  requires  Do  — *  -  ^0  — ^ 

into  (3.10)  finally  gives  that 

^,2^  -  (2  +  =  0, 
dip^ 


.  as  Ip 


The  solution  of  (5.22a)  with  the  required  algebraic  growth  as  tp 
exponential  decay  as  tp  — *  oo  is 
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work  of  tlu'  pn'vious  s.-rtiuii  wh.  io  wr  that  wlu  u  i]  is  in'^ativo  a  similar  critical 

layer  apiJioach  is  reqtiin'il.  see  eijiiatioiis  l  i-lDj. 

T(j  tlescribe  the  detailed  ^tnicnire  we  sul>ln)^^e  that  acrus>,  the  majority  of  tlic  how, 
where  i'  —  0(1),  we  writi' 

1  =  1 1 1  —  ....  f  ( o .  2  o  < ; ) 

and  suppcjse  that  the  streamwi-ie  dejjendt  iice  paramett'r  .i  and  scaletl  i  ^usshow  A  take  the 
for  III 

■)'  --  I  -*  ,  .  .)  am!  )  .  (b.2bi)) 

Sulistitntien  in  i'i.ltj)  yield'^  that  I’o  satmhe'- 

^  -  r„)  -  2A„  (.for  Ao“)  \-,j  =  i-\o.  (5.26) 

This  eiiuatiwii,  which  i-,  identical  to  i  l.iuei.  n.  i -N  ti.  he  ^olve(l  subject  tu  I’o  =  (J  at  r  =  U 
and  I'n  ■  -■  >  o  a^  o  ■  ■  *  oc.  .As  l>.'fuie,  ih' f  -  is  a  c-riri<-al  layer  structure  at  o  =  — .id/Aq  which 
needs  further  stmly  b.  fore  a  full  nuim  rical  sulutidii  of  (.5.26)  >'an  be  eh'ected.  However,  we 
can  draw  the  ceneral  cenclusiun  that  for  <>  "  the  cros.show  nei'ded  to  maintain  neutrally 

stable  Vortices  vaiiishcs  as  v  — 0. 

In  addition  r(j  tie  work  on  neutral  mod«  s  describeil  heri*  wi^  havi'  performed  a  ft'W 
calculations  for  non  neutral  modes.  In  f  iyuii'  (5.S)  we  show  tlu.'  t'h'ect  of  crosshow  on 
the  stabilit}'  of  rim  vortices.  W'e  i-orisidi'red  rlu'  two  freipiem’it'S  =  -1.1  and  took  the 
crosshow  parameter  \  -  .1..2,..3..1  ;ind.  for  various  w.tvi'numbers  k.  solved  (3.10),  (3.11) 
for  the  comj)le.\  valued  .'i .  Fictire  (5.8)  slu-ws  the  cicwth  rate  curves  for  vortices  with  the 
atjove  crosshow  parameters  .  W’e  can  set'  (as  has  already  been  commented  upon)  that  for  a 
hxed  wavi'lencth  vorte.x  an  incr<’asiu”,  crossflow  piocrc-ssively  stabilises  the  how.  Further¬ 
more  we  see  that  instability  occurs  only  for  a  finit*'  ranct'  of  wiliu's  of  the  wavi'iiumber: 
;ind  we  note  thtit  tlm  width  of  this  iinstalde  baiul  <lecreasi'.s  as  the  crossflow  increases.  WT‘ 
also  note  fh.at  tlm  How  is  more  unsttibh'  when  tlm  imiiosed  freqtu'ncy  is  positive'.  WT  now 
conclude  with  seiiim  furtluT  obsei\ations  and  '-(.itim  discussion. 
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§6.  Conclusions. 

In  the  previous  sections  we  have  described  the  process  by  which  the  Gortler  vortex 
instability  mechanism  is  destroyed  by  the  introduction  of  a  w'eak  crossflow  into  a  centrifu- 
gally  unstable  basic  state.  We  choose  to  concentrate  on  the  situation  when  the  crossflow 
is  of  size  Re~  ^  because  that  is  the  crucial  size  at  which  the  crossflow  first  has  a  significant 
effect  on  Gortler  vortices.  In  addition  we  elect  to  concentrate  on  the  large  Gortler  number 
limit  because  that  limit  describes  inviscid  Gortler  vortices  and  nonparallel  effects  can  be 
accounted  for  using  asymptotic  means.  We  found  that  at  order  one  w'avenumbers  the  in¬ 
troduction  of  a  crossflow  causes  the  growth  rate-wavenumber  curve  of  inviscid  theory  to 
split  into  two  distinct  parts.  Thus  above  a  critical  size  of  the  crossflow  a  finite  band  of 
unstable  modes  exists  to  the  right  of  a  =  0.  This  band  of  modes  takes  on  the  asymptotic 
structure  of  crossflow  vortices  when  the  crossflow  is  increased  further.  The  growth  rates 
of  these  modes  increases  monotonically  with  the  crossflow  parameter  and  will  become  for¬ 
mally  of  order  Re^  when  the  crossflow  is  comparable  to  the  streamwise  basic  velocity  field. 
The  second  branch  of  unstable  modes  occurs  at  progressively  higher  wavenumbers  as  the 
crossflow  increases.  The  growth  rates  of  these  modes  at  any  fixed  value  of  the  crossflow 
increases  like  for  large  a  so  it  would  ai)pear  at  first  sight  that,  at  a  fixed  value  of  the 
crossflow,  the  viscous  modes  have  the  highest  growth  rates.  However  this  is  not  the  case 
since  viscous  effects,  as  in  DHS,  serve  so  as  to  stop  this  monotonic  growth  of  f5  with  a.  In 
the  absence  of  any  crossflow  the  growth  rate  peaks  at  wavenumber  of  size  G»  and  has  a 
maximum  of  order  G'»  .  Al  higher  values  of  a  the  growth  rate  tends  to  zero  and  the  right 
hand  branch  neutral  structure  of  Hall  (1982«)  is  retrieved.  In  the  presence  of  a  crossflow 
we  showed  that  the  second,  semi-infinite  band  of  unstable  modes  predicted  on  inviscid 
grounds  becomes  reduced  to  a  finite  band  at  high  wavenumbers  when  the  crossflow  is  for¬ 
mally  of  order  G‘.  At  this  stage  the  growth  rates  of  the  crossflow  vortices  associated  with 
the  modes  to  the  right  of  a  =  0  have  growth  rates  of  size  G»  comparable  with  those  of 
the  viscous  modes.  However  as  the  crossflow  is  increased  further  the  growth  rates  of  the 
crossflow  vortices  increase  whilsi  the  viscous  modes  become  stabilized.  We  further  note 
that  the  unstable  band  of  modes  which  emerge  to  the  right  of  02  connect  with  the  viscous 
modes  in  an  inviscid  manner;  by  this  we  simply  mean  that  their  structure  is  governed  by 
a  critical  layer  scenario  near  02- 

In  Figure  (6.1)  we  have  sketched  the  growth  rate  wavenumber  dependence  determined 
above  for  the  three  situations:  (a)  A*  ~  Gj  .  (b)  A*  ~  G» ,  (c)  A*  G* .  In  (a)  the  crossflow 
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is  strong  enough  only  to  alter  the  inviscid  modes  which  have  growth  rate  of  size  ;  the 
viscous  modes  have  growth  rates  of  order  G«  and  are  not  affected  by  the  crossflow.  In  (b) 
both  modes  have  growth  rates  of  order  G«  but  the  inviscid  modes  are  now  independent  of 
the  wall  curvature  at  zeroth  order;  thus  the  inviscid  modes  hav('  deformed  into  crossflow 
vortices.  The  viscous  modes  at  this  stage  depend  on  the  crossflow  and  become  completely 
stabilized  at  a  finite  value  of  A*/G‘.  However  before  this  occurs  there  are  two  points 
at  which  the  growth  rate  has  a  local  maximum.  The  relative  importance  of  the  maxima 
associated  with  the  viscous  and  inviscid  modes  switches  over  as  the  crossflow  increases.  In 
(c)  only  the  crossflow  vortex  mode  remains  and  at  zeroth  order  it  is  independent  of  the 
curvature  of  the  wall. 

The  question  of  which  of  the  possible  modes  discussed  above  is  the  most  likely  to  be 
observed  in  an  experiment  can  only  be  answered  by  a  receptivity  calculation  of  the  type 
given  by  DHS  for  Gortler  vortices  in  two-dimensional  boundary  layers.  In  fact  the  analysis 
given  in  Secion  6  of  the  latter  paper  is  easily  modified  to  discuss  the  receptivity  problem 
for  the  viscous  modes  described  by  (3.10),  Since  this  type  of  disturbance  is  concentrated 
near  the  wall  then,  following  the  analysis  of  DHS,  it  is  a  relatively  simple  matter  to  show 
that  the  ‘coupling’  coefficient  associated  with  a  roughness  element  at  the  wall  is  0(1).  By 
this  we  simply  mean  that  the  streamwise  velocity  component  of  a  mode  stimulated  by  an 
obstacle  at  the  wall  is  comparable  to  the  fluid  velocity  induced  near  the  obstacle.  How’ever, 
as  in  DHS.  the  coupling  coefficient  depends  on  the  vortex  w’avenumber  and,  if  all  modes 
are  excited,  sufficiently  far  downstream  of  the  obstacle  the  fastest  growing  linear  mode  wdll 
dominate  the  flowfield.  We  sluUl  concentrate  on  the  case  (b)  above  since  by  taking  the 
appropriate  limits  in  that  situation  (a)  and  (c)  can  be  recovered. 

We  suppose  then  that  the  obstacle  stimulates  disturbances  on  a  spanwise  lengthscale 
comparable  to  the  boundary  layer  lengthscale;  this  means  that  the  inviscid  modes  satisfying 
(2.8a)  will  be  stimulated.  We  shall  now  briefly  describe  how  the  coupling  coefficient  for 
this  type  of  mode  can  be  derived. 

Suppose  then  that  the  roughness  element  defined  by 

r  =  /HA',Z),  (6.1) 

generates  a  disturbance  field  governed  by  (2.3).  We  assume  that  there  arc  no  oscilla¬ 
tions  in  the  free  stream  so  that  the  disturbed  flow  will  be  steady.  Following  DHS  it  is 
straightforward  to  show  that  the  system  (2.3)  must  now’  be  solved  subject  to 

U  =  —Ttyh,  V'  =  0,  ir  =  —u'yy.  at  Y  =  0  (6.2) 
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instead  of  =  U'  =  0  at  Y  =  0.  More  precisely  we  note  that  we  must  replace 

ia  by  a  partial  Z  derivative  in  (2.3)  unless  the  X  and  Z  dependences  of  h  are  separable. 
For  siiuphcity  we  asume  that  this  is  the  case  so  that  without  loss  of  generality  we  can 
take  li  to  be  functiiUL  of  X  alone  in  the  hdlowiiig  discussion.  However  we  note  that,  if 
we  wish  to  generate  the  flowfields  in  a  neighbourhood  of  the  obstacle,  the  particular  Z 
dependence  of  h  must  be  accounted  for  when  inverting  the  Fourier  transform  (in  Z)  of  the 
disturbance.  Here  we  are  merely  interested  in  finding  the  coupling  coefficients  associated 
with  stationary  vortices  of  wavelength  comparable  with  boundary  layer  lengthscale  so  it  is 
sufficient  to  take  It  —  h(A’). 

Since  the  forcing  takes  place  at  the  wall  we  must  in  effect  determine  how  an  inviscid 
disturbance  satisfying  (2.8a)  adjusts  in  a  viscous  wall  layer  so  as  to  satisfy  the  correct 
btjundary  con<lirion  at  the  wall.  It  is  easy  to  see  that  this  layer  must  be  of  thickness  G~  « 
so  that  (6.2)  suggest.1  the  following  ex])ansions  for  {U,\\W,P)  in  the  wall  layer  where 

r  =  G^q. 

[U.v. \\\P)  =  (r"  +  .  ,.  +  .  ...G'»P°+  ...)  (l  +  0(G- •))  ,  (6.3) 

The  A'  <lependeuce  of  the  disturbance  fitdd  is  now  taken  care  of  using  a  Laplace  transform, 
suppose  that  h  varies  on  a  G=  huigthscale  in  the  X  direction.  Thus  we  write 

[0.  if  A'  <  A'o, 

'  "  I  /f(Gt[A'  -  A'ol),  if  A  >  Vo. 

and  dett-rmine  the  How  induced  beyond  X  =  A'o-  Thus  we  now  let  3  be  the  Laplace 
transform  variable  with  respt'Ct  to  the  variable  G^ [A"  --  A'o]-  After  a  little  work  we  can 
show  that  in  the  wall  larer  and  are  related  by 

/OO 

Ai{X^q)dq,  (6.4a) 

where  h  is  the  transform  of  //  and 


A  =  3it\i  f  ta/'21- 


wliilst  .\i  is  the  Airy  function.  The  function  1  **  is  then  given  by 

Ai{X^ b)(lb)  <Iq. 


r"  =  -3//.V 


(6.46) 
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w  e  see  that  the  above  solution  has  tendin'^  to  a  constant  when  7  — >  00  so  that  in 
the  main  part  of  the  boundary  layer  the  ap})ropriate  expansions  of  U,  \  IV’  are 

{U,V,W)  =  (G-U'o+  ...,Gh'o+  ...,G'hV’o  4-  ...)  (1  +  0{G-i))  . 


The  function  Vq  is  then  found  to  satisfy  (2.8a):  btit  the  appropriate  boundary  conditions 
for  this  equation  now  become 

I'o  =  ,  at  y  =  0.  Vo - >  0  as  1'  - >  00. 


For  given  values  of  the  transform  variabh's  J  and  a  this  inhomogeneous  differential  system 
for  Vo  can  be  solved.  However  if  d  =  ,3  an  eigenvalue  of  {2.8a.  b.  c)  then  the  solution  of  the 
forced  problem  near  J  must  be  found  by  writing 


V'u  = 


V  bo 

(3  -  3) 


(6.5) 


where  V’oo  =  AIq.  Here  Vq  is  the  eigensolution  of  (2.8a, b.c)  with  3  =  3.  The  constant  A 
may  be  written  as 

Here  the  function  M  is  obtained  by  differentiating  the  h'ft  hand  side  of  ecjuation  (2.8a) 
with  respect  to  3  and  then  setting  V'o  =  Vq.  Having  determined  Vq  we  then  see  that  in 
the  main  part  of  the  boundary  layer  the  disturb(‘d  downstream  velocity  component  in  3,  a 
space  near  3  =  3  given  by 


U 


—  G  ®  A  V  0  Uy 

(;j  —  3  )(.i  a  +  iXaU') 


(6.6) 


It  follows  that  when  the  Fourier-L<aplace  transform  [/(:i,a,  1)  is  inverted  then  the  i)ole  at  3 
will  lead  to  exponentially  growing  eigensoluti.  ,  loportional  to  AG“  ’  (x-Xc)  Thus 
the  coupling  coefficient  for  inviscid  vortex  modes  is  of  order  G~  »  so  that  wall  roughness  is 
a  more  efficient  stimulator  of  the  viscous  modes  which,  from  DHS.  have  an  0(1)  coefficient. 
How’ever  as  the  size  of  the  crossflow  is  increased  the  growth  rates  of  the  inviscid  modes 
increases  to  be  larger  than  the  viscous  gianvth  rates  so  this  slight  diftercuice  in  size  of 
the  couphng  coefficients  will  be  irrelevant  bey<md  a  <’ritical  size  for  the  crossflow  and  the 
observed  instability  will  be  a  Rayleigh  instability  rather  than  a  centrifugal  one. 
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Finally  we  say  a  few  words  about  the  relevance  of  the  above  calculations  to  practical 
flows.  In  the  obvious  practical  situations  where  Gortler  vortices  are  thought  to  be  a  likely 
cause  of  transition  the  basic  state  is  three-dimensional.  Thus  for  example  there  is  little 
doubt  that  the  flow  over  a  turbine  blade  is  three-dimensional  and  our  results  suggest  that 
it  could  therefore  not  support  Gortler  vortices.  On  a  swept  wing  the  basic  state  is  three- 
dimensional  and,  for  order  one  values  of  the  angle  of  sweep,  we  can  be  confident  that  Gortler 
vortices  cannot  exist.  However  though  the  size  of  crossflow  required  to  formally  destroy 
the  mechanism  is  of  size  Re~  ^ G‘  which  is  necessarily  small  because  of  the  definition  of 
the  Gortler  number  the  numerical  constant  multiplying  this  factor  will  be  a  function  of 
the  particular  basic  state.  It  could  well  be  that  for  some  swept-wing  flows  the  Gortler 
mechanism  might  survive  to  modest  angles  of  sweep  appropriate  to  practical  situations. 
Clearly  this  matter  should  be  the  topic  of  a  careful  experimental  investigation  of  the  Gortler 
mechanism  in  three-dimensionu]  flows. 
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igure  (2.1).  The  basic  velocity  profiles  u,  w  defined  by  (2.9 


amplification  parameter  /?  defined  by  (2.6a)  for 


Real  parts  of  the  amplificatioa  parameter  (3  defined  by  (2.6a)  for  various 
crossflow  a)  A  =  0, 1,2,3  and  4,  b)  A  =  4.5,4.75  and  5. 
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Figure  (2.4a,b).  Imaginary  parts  of  the  amplification  parameter  /?  for  various  values  of 
scaled  crossflow  A.  a)  A  =  0,1,2, 3  and  4,  b)  A  =  4.5,4.75  and  5. 


Figure  (2.4r).  Imaginary  parts  ol  the  aiiiplilication  pariaiuefcr  -i  tor  riissiiuw  ^i/• 
A  —  o,  10,  la,  20  ami  25.  F<.>r  i?a<'li  crossflow  tlu-  cross  imiic.atfs  the  classic.il  \.ilui  "t  mu' 
that  (2.12)  IS  satislicd  .ami  so  iifiilr.'illy  stahlc  solutions  ot  the  CjSW  crossIhAs  uist.iliilii 
protili'iii  arise. 
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Figure  (2.5).  a)  Real  parts  and  b)  Imaginary  parts  of  the  amplification  parameter  (3  as 
functions  of  crossflow  A  for  vortex  wavenumbers  a  =  5,  10. 


-0.2 


0.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40  0.45  0.50 

Figure  (4.1).  a)  Real  part  and  b)  Imaginary  part  of /9  as  a  function  of  crossflow  parameter 
A  for  vortices  of  zero  frequency  and  wavenumber  fc  =  1/2. 
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Figure  (4.1).  a)  Real  part  and  b)  Imaginary  part  of  ^  as  a  function  of  crossflow  parameter 
A  for  vortices  of  zero  frequency  and  wavenumber  k  =  1/2. 
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Figure  (5.1).  Values  of  crossflow  A  and  parameter  0  needed  in  o..ier  to  ensure  that 
vortices  of  frequency  fl  =  1  are  neutrally  stable,  a)  A  and  b)  0  as  functions  of  vortex 
wavenumber  k. 
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Figure  (5.2).  Values  of  crossflow  A  and  parameter  P  needed  in  order  to  ensure  that 
vortices  of  frequency  =  0  are  neutrally  stable,  a)  A  and  b)  P  as  functions  of  vortex 
wavenumber  fc. 
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Figure  (5.2).  Values  of  crossflow  A  and 
vortices  of  frequency  A  =  0  are  neutrally 
wavenumber  k. 
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Figure  (5.3).  Values  of  crossflow  A  and  parameter  $  needed  in  order  to  ensure  that 
vortices  of  frequency  n  =  — 1  are  neutraUy  stable,  a)  A  and  b)  0  as  functions  of  vortex 
wavenumber  fc. 
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Figure  (5.4).  Eigenfunctions  U,  V  of  (3.10)  &  (3.11)  for  neutrally  stable  vortices  with 
=  1.35,  n  =  0.  a)  Re(t7),  b)  lm(U),  c)  Re(y),  d)  Iin(F). 
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Figure  (5.4).  Eigenfunctions  U ,  V  of  (3.10)  &  (3.11)  ft 
=  1.35,  n  =  0.  a)  Re(£7),  b)  lm{0),  c)  Re(t^),  d)  Im(y) 
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Figure  (5.5).  Eigenfunctions  U ,  V  of  (3.10)  &  (3.11)  for  neutrally  stable  vortices  with 
=  1.3  X  10-2,  n  =  1.  a)  Re{U),  b)  Iin(C^),  c)  Re(y),  d)  Im(t^). 


66 


,05  0.10  0.15  0.20  0.25  0.30 

Figure  (5.5).  Eigenfunctions  U,  V  of  (3.10)  &  (3.1! 
=  1.3  X  10-2,  n  =  1.  a)  Re(U),  b)  Im(l7),  c)  Re(t^) 
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Figure  (5.6).  Schematic  diagram  of  asymptotic  structure  of  the  solution  of  (3.10)  &c 
(3.11)  for  the  case  of  small  vortex  wavenumber  <C  1.  The  solution  configuration  divides 
into  four  distinct  zones:  Zone  (I)  is  a  thin  region  of  depth  0{k*)  at  a  distance  0{k^)  from 
the  wall,  see  definition  (5.7a).  This  zone  is  embedded  within  (II),  a  region  of  thickness 
0(fc>),  see  (5.12).  The  structure  has  an  additional  thin  viscous  layer  at  the  wall,  Zone 
(///),  which  enables  boundary  conditions  (3.11a)  to  be  met.  Finally,  an  outer  region  (/V) 
of  depth  0(1)  facilitates  the  exponential  decay  of  the  disturbance  solutions  far  from  the 
wall,  see  (5.21). 


67 


1.25  0.30  0.35 

U,  V  of  (3.10)  (3.] 
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Figure  (5.7).  Eigeofuactions  U ,  V  of  (3.10) 
=  8.7  X  10-^  Q  =  -1.  a)  Re(U),  b) 


;ure  (5.7).  Eigenfunctioas  U ,  V  of  (3.10)  &  (3.11)  for  neutrally  stable  vortices  with 
.7  X  10-3,  Q  ^  3^)  Re{f7),  b)  lm(0),  c)  Re(F),  d)  Im(V-). 
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